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This  report  was  prepared  by  ZONA  Technology,  Inc.,  the  prime  contractor,  and  its  Team 
member,  Technosofl,  Inc.  for  AFRLA^ASD,  WPAFB,  Ohio.  It  describes  the  work  performed 
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SECTION  1 


INTRODUCTION 


Summary 

•  Define  Objective  and  Block  Diagram  of  Integrated  Hypersonic  Aerothermoelastic  Program 
Architecture 

•  Background 

•  A  Supersonic-Hypersonic  Vehicle  Design  System:  SHVD 

•  ZONA7U:  A  Unified  Unsteady  Hypersonic/Supersonic  Panel  Method  for  Arbitrary  Wing- 
Body  Configurations 

•  Optimization  Test-Bed  of  ZONA  7U  for  TA  V: ASTROS* 

•  Other  Related  Disciplines  for  TA  V  Design 

•  Description  of  the  following  sections 


This  final  technical  report  describes  work  done  by  the  ZONA  team  (ZONA  Technology,  Inc.  and 
TechnoSofl,  Inc.)  under  SBIR  Phase  I  contract  F33615-01-M-3131,  entitled  integrated 
Hypersonic  Aerothermoelastic  Methodology  for  TAV/TPS  Structural  Design  and  Optimization.  ” 

The  overall  Phase  I  technical  objective  is  to  develop  a  hypersonic  aerothermoelastic 
methodology  for  Trans  Atmospheric  Vehicle  (TAV)/Thermal  Protection  System  (TPS)  structural 
design/optimization  with  a  view  to  integrate  it  with  Adaptive  Modeling  Language  (AML)  (Ref  1) 
into  a  preliminary  TAV  design  software  system. 

Specific  objectives  include: 

•  Establish  interfaces  between  all  key  analysis  software  tools  of  the  preliminary  software 
system  (see  Figure  1.1) 

•  Validate  the  proposed  software  system  by  a  feasibility  study  on  a  selected  TAV 
configuration  (e.g,  X-34) 
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•  Parametric  Geometry 
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FBI  Model 

Mesh  Generator 

1  Aerothermoelas^  Opirrization  \ 

•  Developed  blunt-nose  aerodynamic  methodology,  based  on  a  Strained-Coordinate 
technique  analytically  applied  to  the  local  panel  by  matching  the  Chernyi’s  similarity 
solution  (Ref  2),  including  Lees’  hemisphere  solution  (Ref  3),  with  ZONA’s  unified 
hypersonic/supersonic  pulsating-cone  solution  (Ref  4)  on  downstream  panels 

•  Seamlessly  integrated  the  SHARP  module  of  the  MARKV  code  (Ref  5)  into  ZONA7U 
for  aeroheating  analysis 

•  Developed  a  finite  element  based  streamline  code  called  ZSTREAM  that  adopts  the 
inviscid  surface  velocities  generated  by  ZONA7U  as  input  to  yield  high  quality 
streamline  solutions 

•  Integrated  the  EXITS  module  of  MINIVER  (Ref  6)  with  ZONA7U+SHABP  for  TPS 
sizing 

•  Demonstrated  the  integrability  and  the  trim  analysis  capability  of  ASTROS*  for  a 
flexible  X-34  in  hypersonic  maneuver  and  re-entry  phase. 

•  Validated  Central  Methodologies  (Blocks  1-5  in  Figure  1.1)  required  for 
aerothermodynamic  optimization  individually  for  X-34,  a  selected  TAV  configuration  for 
methodology  demonstration. 
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1.1  Bacl^round 

Aerothermoelastic  analysis  has  become  a  required  discipline  for  TAV  design. 
Aerothermoelasticity  is  a  synergic  disciplinary  of  aerothermodynamics  and  aeroelasticity.  In  a 
hypersonic  extreme  environment,  aerothermoelastic  effects  will  strongly  influence  the  TPS 
sizing  and  the  integrated  TAV/TPS  structural  design.  During  the  hypersonic  flight  phase,  the 
aerothermoelastic  load  will  cause  TAV  deformation,  which  in  turn  will  impact  the  structural 
integrity  of  the  TAV/TPS  system.  If  designed  properly  with  the  main  structures,  the  TPS  will 
serve  as  a  part  of  the  load-carrying  structure,  thus  helping  further  reduce  the  total  weight  of  the 
TAV.  For  this  reason,  both  the  TAV  and  its  TPS  require  an  accurate  aerothermodynamic  loads 
prediction  method  to  couple  with  an  optimization  method  in  order  to  achieve  a  viable  TAV/TPS 
structural  design. 

Other  important  disciplines  that  will  influence  TAV/TPS  design  are  the  trajectory  analysis,  the 
TPS  sizing  and  thermal  analysis.  To  integrate  all  these  disciplines  and  turn  them  into  a 
Multidisciplinary  Design  and  Optimization  (MDO)  design  tool  for  TAV/TPS  presents  great 
challenges.  A  series  of  Computational  Fluid  Dynamics  (CFD)  works  and  test  data  were  directed 
towards  the  X-34  research  and  developed  by  NASA  and  Orbital  Sciences  Corporation  in  the  last 
decade  (Refs  7-12).  Adopting  the  CFD  approach  for  an  effective  TAV/TPS  design  tool  in  a 
hypersonic  environment  is  prohibited  by  the  slow  CFD  turn  around  time.  On  the  other  hand, 
other  more  expedient  computational  methods  (e.g..  Refs  13-14)  utilized  in  treating  all  these 
required  disciplines  have  been  developed  individually  to  a  certain  extent  including  the  focused 
aerothermodynamic  program.  Nonetheless,  all  these  previous  computational  approaches  lack  a 
main  design-oriented  program  with  data  management  capability  and  multidisciplinary 
design/optimization  perspective.  An  MDO  oriented  program  such  as  ASTROS*  (Automated 
STRuctural  Optimization  System)  or  NASTRAN  should  be  the  central  piece  of  a  valid  TAV/TPS 
design  tool.  If  ASTROS*  is  selected  then  a  compatible  aerothermodynamic  program  must  be  a 
high-fidelity  one  in  order  to  interface  with  a  structural  FEM  module.  This  requirement  will 
probably  rule  out  the  existing  efficient  but  non-FEM  compatible  types  of  aerodynamic  prediction 
programs  such  as  APAS  (Ref  15),  Datcom  (Ref  16),  or  AP98  (Ref  17).  Clearly,  a  high  fidelity, 
computationally  efficient  hypersonic  aerothermoelastic  methodology  is  lacking. 

1.2  A  Supersonic-Hypersonic  Vehicle  Design  System;  SHVD  (Ref  18) 

Adaptive  Modeling  Language  (AML),  developed  by  TechnoSoft,  Inc.,  offers  the  advanced, 
object-oriented  engineering  modeling  language  to  enable  the  modeling  and  simulation  of  the 
entire  product  development.  Based  on  AML,  TechnoSoft  and  Lockheed  Martin  are  developing  a 
Supersonic-Hypersonic  Vehicle  Design  (SHVD)  system,  which  is  an  object-oriented,  web- 
enabled  distributed  fi-amework  environment  for  design  analysis  and  simulation  of  TAV.  The 
SHVD  system  automates  and  manages  the  data  transfer  between  various  design,  analyses  and 
simulation  tools,  including  aerodynamics,  aero-heating.  Thermal  Protection  System  (TPS), 
propulsion  system,  trajectory  analysis,  structural  weight  optimization  and  cost.  The  SHVD 
system  builds  upon,  and  leverages  years  of  software  development  and  TAV  design  domain 
knowledge  offered  by  the  entire  SHVD  team:  1)  AML,  a  web-enabled  Adaptive  Modeling 
Language  from  TectmoSoft;  2)  IMD,  an  Interactive  Missile  Design  system  from  Lockheed 
Martin  Missile  and  Fire  Control;  3)  Numerous  aerodynamic,  aero-thermal,  propulsion  and 
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trajectory  codes  and  detailed  domain  knowledge  offered  by  Lockheed  Martin  Aeronautics 
Company  (Forth  Worth,  Skunk  Works  and  Marietta);  4)  Structural/TPS  sizing,  optimization 
applications  and  domain  knowledge  offered  by  Collier  Research;  5)  Software  support/validation 
for  TAV  design  by  NASA/Langley  and  NASAMarshall. 

The  SHVD  development  resulted  in  an  advanced  multidisciplinary  capability  for  TAV 
design/simulation  under  hypersonic  extreme  environment;  the  enabling  software  methodologies 
include; 

-  Aerodynamic  Analysis:  S/HABP,  APAS,  PANAIR,  VUAERO 

-  Thermal  Analysis:  MINIVER,  FEM/SINDA 

-  Boundary  Layer/ Aero-Heating  Analysis:  S/HABP  MarkV 

-  Trajectory  Analysis:  POST 

-  Structural  Analysis:  NASTRAN 

For  realistic  TAV  design/analysis  in  an  extreme  hypersonic  environment,  aero-heating  problems 
related  to  aerothermodynamic  and  aerothermoelasticity  will  strongly  impact  the  TAV/TPS 
structural  design  compatibility,  hence  the  vehicle  structural  integrity.  Further,  the  aeroelastic 
instability  induced  by  aero-heating  could  lead  to  serious  divergence  or  flutter  problems  in  its 
hypersonic  flight  phase.  Close  examination  of  the  aerodynamic  software  capability  of  SHVD 
(e.g..  Ref  1 8)  reveals  that: 

-  S/HABP-APAS  does  not  have  unsteady  aerodynamic  capability  needed  for  flutter, 
divergence  and  ASE  instability  analysis 

-  It  does  not  generate  Aerodynamic  Influence  Coefficient  (AIC)  matrices  needed  for 
sensitivity  analysis  to  avoid  repetitive  aerodynamic  computation  in  an  optimization 
procedure 

-  It  is  not  a  PANAIR-level  high-fidelity  panel  method;  hence,  it  offers  no  solution 
refinement  from  conceptual  to  preliminary  design  stage 

On  the  other  hand,  unsteady  aerodynamics-structural  dynamics  coupling  via  Panel-FEM 
interface  has  been  a  widely  practiced  methodology  in  both  loosely  or  tightly-coupled  levels. 
Only  until  the  recent  ZONA  unsteady-hypersonic  development,  such  a  Panel-FEM  methodology 
for  hypersonic  TAV  design  has  not  been  available  due  to  the  lack  of  a  suitable  hypersonic 
aerodynamic  panel  method. 

For  decades,  the  difficulty  that  has  hampered  the  development  of  a  viable  hypersonic  panel 
method  lies  in  the  following: 

-  Superinclined  Panel:  Supersonic  kernel  integral  becomes  singular  at  a  Mach  number  at 
which  the  inclined  panel  slope  exceeds  that  of  the  Mach  wave.  Computation  breaks  down 
beyond  this  Mach  number 

-  Flow  Rotationality:  Strong  shock  induced  flow  rotationality  carmot  be  included  in  the 
conventional  potential  flow-based  panel  method.  Accuracy  will  deteriorate  rapidly  as  the 
supersonic  Mach  number  increases  toward  the  hypersonic  range 
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1.3  Z0NA7U:  A  Unified  Unsteady  Hypersonic/Supersonic  Panel  Method  for 
Arbitrary  Wing-Body  Configurations 

With  continuous  R&D  in  hypersonics  since  1995,  ZONA  Technology,  Inc.  (ZONA)  has  made 
major  breakthroughs  in  overcoming  the  above  two  issues.  The  result  is  a  unique  software  product 
ZONA7U,  a  high-fidelity  unified  unsteady  hypersonic  panel  method  (Ref  4).  ZONA7U  (U 
stands  for  unified)  has  the  following  capabilities: 

-  It  is  a  fi-equency-domain  and  s-domain  aerodynamic  for  unified  supersonic/hypersonic 
Mach  numbers  up  to  the  Newtonian  limit,  hence  capable  of  periforming  hypersonic 
aeroelasticity/aerothermoelastic  applications  including  flutter,  divergence,  gust  and 
aeroservoelasticity  (ASE)  instability 

-  It  provides  unified  supersonic/hypersonic  AIC  matrix,  hence  is  readily  applicable  for 
structural  optimization  procedure  with  aeroelastic  constraints 

-  It  can  aerodynamically  model  complex  air  vehicles  such  as  conventional  aircrafts, 
blended  wing-body  and  TAV  configurations,  hence  a  high-fidelity  panel  method 

Based  on  the  formulation  of  strained  coordinates  in  conjunction  with  the  local  pulsating  body 
analogy,  ZONA7U  can  accurately  approximate  the  nonlinear  thickness  effect  and  the  shock- 
induced  rotationality  in  the  unified  supersonic/hypersonic  Mach  range  up  to  the  Newtonian  limit. 
Both  of  these  effects  are  ignored  by  the  linear  theory  and  are  overestimated  by  the  Piston  theory 
(Refs  19,  20,  21).  In  the  following  examples,  it  can  be  seen  that  the  ZONA7U  results  agree  well 
with  the  Euler  solutions  for  various  classes  of  wings  and  bodies. 

•  ZONA7U for  Various  Wing Planforms 


Figures  1.2-1. 5  present  the  ZONA7U  unsteady  damping  derivative,  generalized  aerodynamic  forces  and 
flutter  solutions  for  various  wing  planforms  at  hypersonic  Mach  numbers  showing  comparable  accuracy. 


(a)  (b)  (c) 

Figure  1.2  ZONA7U  Damping-in-Pitch  Derivatives  of  a  Rectangular  Wing  with  a  Diamond 
Profile  versus  Airfoil  Thickness  (o)  at:  (a)  M=2.0,  (b)  M=5.0,  and  (c)  M=10.0. 
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Figure  1.3  Z0NA7U  Damping-in-Pitch  Derivatives  of  a  Rectangular  Wing  with  a  Diamond 
Profile  versus  Mach  Number  (cf=15°,  h/c=0.5). 


Figure  1.4  Effect  of  Reduced  Frequency  on  Generalized  Aerodynamic  Forces  for  Oscillating 

Panels  ZONA7U  at  M  =  5.0,  a  =  2°,  N  =  2 
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•  Z0NA7U for  Bodies 

Figures  1.6-1.15  present  the  Z0NA7U paneling,  pressure  distributions  and  unsteady  stability  derivatives 
for  various  bodies  showing  excellent  agreement  with  CFLSD/Euler  solutions  (Ref  22)  and  measure  data. 


Under  a  recent  Army/REDC  support, 
ZONA  has  further  extended  ZONA7U  to 
treat  body-fin  configurations  at  Mach  6.0. 
To  circumvent  the  superinclined  panel 
problem  (i.e.,  when  the  Mach  line  cuts  into 
the  body  panel  due  to  high  Mach  number), 
we  introduce  an  equivalent  Mach  number 
transformation  to  recast  the  physical 
problem  into  a  new  coordinate,  whereby 
the  body  undergoes  a  compressibility 
stretch  in  the  axial  direction. 


(a)  ZONA7U Panel  Model  (b)CFL3DMesh 

Figure  1.6  ZONA7U  and  CFL3D  Models  of  the 
CKEMBody 


For  the  inclusion  of  flow  rotationality  effects,  we  have  established  a  local  pulsating  body  analogy 
to  extend  Sims’  (Ref  23),  Brong’s  (TRef  24),  and  Dorodnitsyn’s  (Ref  25)  steady  Euler  solutions. 
A  detailed  theoretical  formulation  of  ZONA7U  can  be  found  in  References  4  and  26  (Chen  and 
Liu). 


Note  that  CFL3D  requires  over  2  hours  of  computer  time  for  each  bent-nose  case  whereas 
ZONA7U  takes  only  1  minute. 


Undeflected  Case,  Body  &  Nose  at  the  same  AOA 
(Nose  -  Body  =  -2®,  atM-6.0) 
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Figure  1.7  ZONA7U  Pressure  Distributions  and  Aerodynamic  Force/Moments  along  the  CKEM 
Body  at  M  -  6.0  for  Various  Bent-Nose  Angles  and  Angles  of  Attack 
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Cp/CpO 


Figure  1.9  Pressure  Distributions  on  a  Henrisphere  Body  at  M  =  5.8 
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Figure  1.12  Variation  with  Mach  Number 
of  the  Stability  Derivative  (at  Cone  Angle  = 
10») 


Figure  1.13  Variation  with  Mach  Number 
of  the  Stability  Derivative  (at  Cone  Angle  = 
20°) 
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Figure  1.14  ZONA7U  Damping-in-Pitch  Figure  1.15  ZONA7U  Damping-in-Pitch 
Moment  Coefficients  for  a  Cone  Frustum  Moment  Coefficients  for  an  Ogive  Cylinder  at 

Various  Pivot  Locations 


1.4  Optimization  Test-Bed  of  ZONA7U  for  TAV:  ASTROS* 

Under  recent  Army/RDEC  support,  ZONA  has  successfully  integrated  ZONA7U  with 
ASTROS*  and  performed  structural  optimization  for  a  hypervelocity  Compact  Kinetic  Energy 
Missile  (CKEM)  design  with  flexible  body  control  at  Mach  6.0  (Ref  26).  ASTROS*  is  ZONA 
Technology’s  enhanced  version  of  ASTROS  that  is  seamlessly  integrated  with  ZONA’s  Unified 
Aerodynamic  ZEARO  system  and  Aeroservoelasticity  (ASE)  module.  ASTROS  (Automated 
STRuctural  Optimization  System)  is  a  proven  engineering  design/analysis  and  optimization 
software  which  includes  most  of  the  essential  aerospace  disciplines  that  impact  an 
aircraft/missile  structural  design  (isotropic  or  composites). 

Under  a  two-year  contract  support  by  AFRL,  ZONA  Technology  Inc.  (ZONA)  has  further 
enhanced  the  software  system  by  seamlessly  integrating  several  additional  engineering  modules 
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into  ASTROS  (Refs  27,  28).  Figure  1.16  shows  all  the  essential  modules  of  ASTROS*, 
including  die  ZAERO  aerodynamic  module,  the  smart  structure  module,  the  ASE  module,  the 
trim  module  and  the  aeroelastic  stability  module.  The  functionalities  and  features  of  each 


module  are  shown  in  Figure  1.16.  Note  that  ASTROS*  =  ASTROS  +  ZAERO  (Refs  29, 30) 


Smart  Structures  Module 


Modeling  of  PZT  and  SMA 
activations  and  computes 
the  Induced  aerodynamic 
control  forces. 


ZAERO  Aerodynamic 
Module 


Unified  steady/unsteady 
aerodynamics  for  subsonic, 
transonic,  supersonic,  and 
hypersonic  flows. 


Sensitivity  Analysis 


Provides  sensiti\^ty  of 
stresses,  stability,  and 
perfonnance  with 
respect  to  structural 
design  variables. 


Structural  Finite  Element  Module 


NASTRAN-compatible  FEM  analysis 


Trim  module 


Trim  Module 


Static  aeroelastic 
analysis  for  flight 
loads  and  optimum 
trim  solutions 


Aeroelastic  Stability 
Module 


Provides  true  damping 
flutter  solutions. 


Optimization  Module 

\  An  optimizer  driving  all 
modules  to  achieve  the 
optimum  design. 


ASE  Module 


State-space  aeroservoeiastic  analysis  including  S  ISO/MI  MO 
control  system  for  stability  and  gust  loads  analysis. 


Figure  1.16  Engineering  Modules  in  ASTROS* 


ZONA  plans  to  use  ASTROS*  as  the  ZONA7U  test-bed  for  TAV  structural  design/optimization. 
Once  shown  successful,  ZONA7U  compatibility  is  assured  for  subsequent  interface  with  other 
structural  FEM  codes  such  as  NASTRAN. 


1.5  Other  Related  Disciplines  for  TAV  Design 

Aerothermoelastic  analysis  is  a  major  concern  for  load-carry  structures  and  Thermal  Protection 
System  (TPS)  structural  design  of  TAV.  In  addition  to  the  strength  and  flutter/divergence 
constraints,  the  thermal-stress  and  the  material  property  degradations  due  to  the  aero-heating 
effects  must  also  be  taken  into  an  account  for  designing  die  load-carry  structures. 

The  TPS  sizing  is  normally  dominated  by  the  boundary  conditions  of  the  thermal  analysis. 
Petley  et  al.  (Ref  31)  outlined  an  excellent  design  procedure  of  TPS  sizing  in  which  they 
conducted  a  trajectory  analysis  to  obtain  a  transient  heating  profile  over  a  Mach  10  cruise 
vehicle.  Two  types  of  TPS  design  concepts  were  presented;  one  for  the  vehicle  skin  over  the 
integral  tank  structure  and  the  other  for  skin  with  no  integral  tank  to  provide  the  hydrogen  heat 
sink  capability.  For  TPS  sizing,  SHVD  has  MINIVER  and  FEM/SINDA  as  software  tools  for 
design/analysis. 
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Trajectory  simulation/optimization  is  required  for  aero-heating  analysis  and  vehicle  mission 
closure  for  TAV.  POST  (Program  to  Optimized  Simulated  Trajectories)  developed  by  NASA  is 
included  in  SHVD  for  fuel  minimization  and  other  objectives  for  a  given  mission. 
Aerothermoelastic  analysis  requires  local  temperature  gradient  heat  transfer  and  heat  rate  over 
the  outer  mold-line  (OML)  of  the  TAV.  This  in  turn  requires  the  hypersonic  inviscid  and 
compressible  boundary  layer  interaction.  ZONA7U  will  provide  the  edge  pressure  and  velocity 
which  will  be  coupled  with  the  boundary  layer/aero-heating  module  of  S/HABP  MarkV  to 
replace  MarkV’ s  inviscid  module.  The  output  is  the  above  local  heat  transfer  quantities,  heat 
transfer  coefficients  and  equivalent  skin  friction  coefficients.  The  last  parameter  will  be  a  new 
trajectory  input  from  the  integrated  ZONA7U. 

1.6  Description  of  the  Sections 

There  are  9  sections  in  the  present  report. 

•  Section  2  describes  our  methodology  development  of  the  hypersonic  aerodynamics  and 

aerothermodynamics  for  basic  vehicle  component. 

•  Section  3  describes  our  development  of  a  surface  streamline  routine,  ZSTREAM,  which  is 
essential  in  aerodynamics/aeroheating  interface,  and  an  integral  part  of  the  present 
methodology. 

.  Section  4  describes  the  application  of  the  present  hypersonic  aerothermodynamics 
methodology  to  X-34,  a  selected  TAV  configuration. 

.  Section  5  describes  how  we  employ  ASTROS*  to  perform  analysis  for  a  flexible  X-34, 
under  aerothermodynamic  loading,  demonstrating  the  integrability  of  the  present 
methodology  with  ASTROS*  as  a  design/optimization  central  software. 

•  Section  6  describes  the  elementary  TPS  sizing  procedure  in  anticipation  to  a  proposed 
automated  optimized  procedure  (See  Section  8). 

•  Section  7  describes  aeroheating  history  of  the  stagnation  point  of  X-34  throughout  two 
different  trajectories  and  shows  validation  with  that  due  to  MINIVER. 

•  Section  8  is  the  phase  I  concluding  researches/recommendation. 

•  Section  9  describes  the  proposed  phase  II  plan. 

The  interfacing/interaction  plan  of  AML  without  present  hypersonic  aerothermodynamics 
methodology  is  shown  in  the  Appendix  A. 
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SECTION  2 


HYPERSONIC  AEROTHERMODYNAMICS  DEVELOPMENT: 
POINTED  AND  BLUNT  BODIES 


Summary 

In  this  section,  we  describe  ZONA ’s  extended  development  of  its  unified  hypersonic/supersonic 
aerodynamic  tool  ZONA  7U  to  aerothermodynamic  applications.  In  so  doing,  a  new  surface 
streamline  method  has  been  developed  based  on  ZONA  7U.  ZSTREAM,  the  ZONA  developed 
streamline  scheme,  is  Mach  number  dependent  and  uniformly  valid  everywhere  including  the 
stagnation  point.  ZSTREAM  is  to  replace  the  role  of  QUADSTREAM  in  SHARP.  Thus  ZONA 
unified  hypersonic/supersonic  aerothermodynamic  methodology  is  developed  by  suitably 
coupling  ZONA  7U  +  ZSTREAM  +  SHARP.  Two  cases  studied  in  hypersonic 
aerothermodynamics  for  pointed  body  (CKEM)  and  for  a  15°  blunt  cone  are  shown  to 
demonstrate  the  ZONA  developed  methodology.  Our  computed  results  in  Cp  and  heat  flux  are 
validated  with  those  obtained  by  a  CFD  method  (CFL3D  +  LATCH,  also  worked  out  by  ZONA). 
Overall  good  correlations  are  found  except  that  the  computing  efficiency  of  the  ZONA  method  is 
one  to  two  orders  faster  than  that  of  the  CFD  method. 


2.1  ZONA7U  Unified  Hypersonic  Aerodynamics  (Block  1) 

Previous  description  of  ZONA7U  capability  includes: 

-  It  generates  unsteady  aerodynamics,  aeroelasticity  and  aeroservoelasticity  instability 
solution  for  design  constraints 

-  It  provides  Unified  Mach  number  AIC  Matrix  readily  to  be  integrated  with  structural  FEM 
of  ASTROS*  for  MDO 

-  It  is  generally  applicable  to  arbitrary  wing-body  configuration 

In  addition  to  the  above,  ZONA7U  can  handle  sharp  nose  as  well  as  blunt  nose  bodies  at 
arbitrary  angles  of  attack  up  to  flow  separation.  The  blunt-nose  aerodynamic  methodology  is 
based  on  a  Strained-Coordinate  technique  analytically  applied  to  the  local  panel  by  matching  the 
Chemyi’s  similarity  solution,  including  Lees  hemisphere  solution,  with  ZONA’s  unified 
hypersonic/supersonic  pulsating-cone  solution  on  downstream  panels.  In  Phase  I,  a  validation 
effort  of  the  ZONA7U  pressure  solution  was  extensive  showing  good  correlation  with  test  data 
and  CFD  solution  up  to  flow  separation  (see  Figs  2. 1  -  2. 13). 

2.2  ZSTREAM  Development  for  Aerothermodynamics  (Block  2) 

In  Phase  I,  we  have  seamlessly  integrated  the  SHABP  module  of  the  MARKV  code  into 
ZONA7U  for  aeroheating  analysis. 

SHABP  performs  the  boundary  layer  /  aeroheating  computations  along  each  streamline  using  the 
empirical  equations  documented  in  the  first  Phase  I  progress  report.  In  the  original  SHABP 
module,  the  Newtonian  steepest  decent  method  (called  “QUADSTREAM”)  is  employed  for 
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streamline  computation.  However,  it  was  found  that  the  Newtonian  steepest  decent  method 
experienced  numerical  problems  on  complex  configurations  that  contain  highly  twisted  panels. 
Because  the  four  comer  points  of  a  twisted  quadrilateral  panel  are  not  located  on  the  same  plane, 
a  gap  exists  between  two  adjacent  twisted  panels.  The  streamline  tracing  procedure  of  the 
QUADSTREAM  method  is  not  robust  enough  to  overcome  this  problem. 

Another  shortcoming  of  QUADSTREAM  is  its  independence  of  Mach  number.  The  only 
information  needed  to  generate  the  streamline  is  the  freestream  velocity  vector  and  a 
quadrilateral  discretization  of  the  geometry.  Obviously,  the  accurate  streamline  locations  should 
be  determined  by  the  local  velocity  vectors  over  the  surface  of  the  configuration.  Because  the 
aerodynamic  methods  in  the  SHABP  code  are  all  based  on  empirical  methods,  they  do  not 
generate  the  local  velocity  vectors  over  the  surface;  only  the  local  pressure  is  computed. 
However,  this  is  not  the  case  if  the  ZONA7U  code  is  used  because  of  its  hypersonic  panel 
method  formulation  that  generates  the  local  velocity  vectors. 

In  Phase  I,  we  have  developed  a  finite  element  based  streamline  code  called  ZSTREAM  that 
adopts  the  inviscid  surface  velocities  generated  by  ZONA7U  as  input  to  yield  high  qu^ity 
streamline  solutions.  The  details  of  the  theoretical  formulation  of  ZSTREAM  will  be  described 
in  Section  3  in  detail. 

To  validate  the  ZONA7U/ZSTREAM/SHABP  procedure,  we  have  performed  the  aeroheating 
analysis  on  three  configurations,  namely  the  CKEM  (Compact  Kinetic  Energy  Missile)  body  at 
M=6.0  and  a=2°,  a  15°  blunt  cone  at  M=10.6  and  a=0°,  5°,  and  10°  and  the  X-34  wing-body 
configuration  at  M=6.0,  a=9°  and  15.22°.  These  results  are  compared  to  the  CFD  results  using 
CFL3D  +  LATCH  (Ref  25).  These  comparisons  on  the  inviscid  Cp  and  heat  transfer  rates  are 
shown  as  follows: 

2.3  Case  Study  (A);  CKEM  Body  at  M  =  6.0  and  a  =  2° 

The  CKEM  body  consists  of  a  sharp  ogive  nose  and  a  cylinder  body.  Fig  2. 1  shows  the  inviscid 
surface  pressure  distribution  computed  by  ZONA7U  and  the  CFL3D  Euler  solver.  Fig  2.2 
depicts  the  comparison  of  the  wind-side  Cp  distribution  computed  by  ZONA7U  and  CFL3D.  In 
both  figures,  excellent  agreement  between  the  ZONA7U  results  and  the  CFL3D  results  can 
clearly  be  seen.  The  streamlines  computed  by  ZSTREAM  using  the  ZONA7U  generated  surface 
velocities  are  shown  in  Fig  2.3 .  For  clarity,  only  streamlines  associated  with  panels  at  the  rear  of 
the  CKEM  body  are  shown.  Based  on  these  ZSTREAM  results  and  the  ZONA7U  inviscid  Cp 
results,  the  Laminar  heat  transfer  rates  (q)  distribution  is  shown  in  Fig  2.4  and  compared  with 
results  computed  by  the  CFL3D  +  LATCH  code.  It  should  be  noted  that  the  streamline 
computation  procedure  of  LATCH  is  based  on  an  integral  method  that  contains  a  singularity  at 
the  stagnation  point.  Due  to  this  singularity,  LATCH  cannot  provide  the  heat  transfer  rate  near 
the  stagnation  point  resulting  in  a  “cut-out”  of  q  at  the  nose.  The  comparison  of  q  on  the  wind- 
side  surface  between  ZONA7U  +  SHABP  and  CFL3D  +  LATCH  using  the  boundary  condition 
Tw  =  540°R  is  shown  in  Fig  2.5.  At  the  ogive  nose  region,  ZONA7U  +  SHABP  slightly  over¬ 
predicts  the  heat  transfer  rates  when  compared  to  CFL3D  +  LATCH.  This  discrepancy  is 
probably  caused  by  the  shortcoming  of  the  empirical-equation-based  methodology  in  SHABP. 
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ZONA71!  CFL3D/Euler 


Figure  2.1  Inviscid  Surface  Pressure  Distribution, 

Sharp  Cone/CKEM:  M,  =6.0,a=2°,  p, -2.66  lb/ft%  T,  =89.97rR,  =540<’R 


Figure  2.2  Wind-Side  Inviscid  Surface  Pressure  Distribution  ((1)=180®). 
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Figure  2.3  Streamlines  Computed  by  ZX>NA7U/ZSTREAM,  M=6.0. 


Z»NA7U  +  SHABP 


CFUD/Euler + LATCH 


Figure  2.4  Laminar  Heat  Transfer  Rates  (Btu/ft  -s). 
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Figure  2.5  Wind-Side  Laminar  Heat  Transfer  Rates  ((j)-180®). 


2.4  Case  Study  (B)  15®  Blunt  Cone  at  M=10.6,  a=0°,  5°,  and  10® 

The  15°  blunt  cone  geometry  consists  of  a  15°  cone  with  a  round  nose  of  1.1”  radius.  The 
inviscid  surface  Cp  distributions  computed  by  ZONA7U  and  CFL3D  are  presented  in  Fig  2.6  for 
a  =  0°,  2.7  for  a  =  5°,  and  2.8  for  a  =  10°  where  good  correlations  near  the  nose  region  can  be 
seen.  Note  that  at  a  =  0°,  the  CFL3D  result  shows  some  numerical  oscillation  while  ZONA7U 
remains  smooth.  The  comparisons  of  the  inviscid  Cp  between  the  test  data  and  the  computed 
results  of  ZONA7U  and  CFL3D  are  depicted  in  Figs  2.9,  2.10,  and  2.11  for  the  a  =  0°,  5°,  and 
10°  cases,  respectively.  At  a  =  0°,  both  the  ZONA7U  and  CFL3D  results  correlate  well  with  the 
test  data.  At  a  =  5°  and  10°,  the  CFL3D  results  give  a  better  correlation  with  the  test  data  than 
ZONA7U.  This  is  expected  since  CFL3D  solves  the  Euler’s  equations;  whereas,  ZONA7U 
suffers  from  the  attached-flow  assumption  and  consequently  loses  its  accuracy  at  high  angles-of- 
attack.  However,  CFL3D  requires  much  longer  computing  time  than  ZONA7U.  For  the  present 
case,  CFL3D  requires  30  hours  of  computing  time  whereas  ZONA7U  requires  only  10  minutes. 

Based  on  the  surface  velocities  generated  by  ZONA7U,  the  streamline  computed  by  ZSTREAM 
at  a  =  0°,  5°,  and  10°  are  presented  in  the  Fig  2. 12.  The  Laminar  heat  transfer  rates  computed  by 
ZONA7U  +'SHABP  and  CFL3D  +  LATCH  at  a  =  0°,  5°,  and  10°  are  depicted  in  Figs  2.13, 
2.14,  and  2.15,  respectively.  Again,  due  to  the  singularity  in  the  integral  method,  the  CFL3D  + 
LATCH  cannot  provide  q  at  the  stagnation  point  near  the  nose  whereas  ZONA7U  +  SHABP 
does  not  have  such  a  problem.  Figs  2.16,  2.17,  and  2.18  present  the  comparison  of  q  between 
the  computed  results  of  ZONA7U  +  SHABP,  of  CFL3D  +  LATCH  and  the  test  data  for  the  a  = 
0°,  5°,  and  10°  cases,  respectively.  The  temperature  on  the  wall  is  fixed  at  540°R.  In  general, 
the  CFL3D  +  LATCH  results  correlate  better  with  the  test  data  than  those  of  the  ZONA7U  + 
SHABP  results.  Again,  the  discrepancy  of  the  ZONA7U  +  SHABP  is  probably  due  to  the 
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inaccurate  Cp  at  high  angles  of  attack  and  the  shortcoming  of  the  empirical-equation-based 
methodology  in  SHABP. 


a  =  0° 

ZONA7U  CFL3D/Euler 


F^re  2.6  Inviscid  Surface  Pressure  Distribution 
IS**  Blunt  Cone  at  M„  =10.6,  a=0°,  p„=2.66  Ib/ft*,  T„=89.97rR,  T^=540“R 
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a  =  5« 

ZONA7U  CFL3D/Euler 


Figure  2.7  Inviscid  Surface  Pressure  Distribution 


15®  Blunt  Cone  at  M„=10.6,  a=5°,  p„=2.66  lb/ft^  T„=89.971®R,  =540®R 


a  =  10® 

ZONA7U  CFL3D/Euler 


Figure  2.8  Inviscid  Surface  Pressure  Distribution 


15®  Blunt  Cone  Case  at  M„  =10.6,  a=10®,  p«=2.66  Ib/fP,  T,=89.971®R,  T,  =540®R 
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Figure  2.10  Inviscid  Surface  Pressure  Distribution  (<|)=180°) 
on  a  15°  Blunt  Cone  at  M=10.6  and  a=5°. 
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Figure  2.1 1  Inviscid  Surface  Pressure  Distribution  ((t)=180®) 
on  a  15®  Blunt  Cone  at  M=10.6  and  a=10®. 


Figure  2.12  Streamlines  on  a  15  Blunt  Cone  at  M=10.6;  (a)  a=0°,  (b)  a=5°,  (c)  a-10°. 
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a  =  0° 


CFL3D/Euler  + LATCH 


ZONA7U  +  SHABP 


Figure  2.13  Laminar  Heat  Transfer  Rates(Btu/fP-s) 

15«  Blunt  Cone  at  M„=10.6,  a=0°,  p„=2.66  Ib/fP,  T„=89.97rR,  T,  =540»R 


a  =  5® 

ZONA7U  +  SHABP  CFL3D/Euler  +  LATCH 


15“  Blunt  Cone  at  =10.6,  a=5°,  p„=2.66  Ib/tf ,  T„=89.971“R,  =540“R 
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a  =  10° 

Z0NA7U  +  SHABP  CFL3D/Euler  +  LATCH 


Figure  2.15  Laminar  Heat  Transfer  Rates(Btu/fP-s) 

15”  Blunt  Cone  at  M„  =10.6,  a=10“,  p„=2.66  Ib/ft^,  T„=89.971”R,  T„  =540”R 


Figure  2.16  Laminar  Heat  Transfer  Rates  on  a 
15”  Blunt  Cone  at  M„  =10.6,  a=0°,  p„=2.66  Ib/ft^,  T„=89.971”R,  =540”R 


Figure  2.17  Laminar  Heat  Transfer  Rates  on  a 
15“  Blunt  Cone  at  M„  =10.6,  a=5°,  p„=2.66  Ib/ft",  T„=89.971“R,  T„=540“R 


Figure  2.18  Laminar  Heat  Transfer  Rates  on  a 
15“  Blunt  Cone  at  M„=10.6,  a=10°,  p„=2.66  rf,  T„=89.971“R,  T„=540“R 
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SECTION  3 


ZSTREAM  for  Hypersonic  Aerothermodynamic  Methodology 


Summary 

The  ZONA  development  of  ZSTREAM  was  prompted  by  the  breakdown  of  QUADSTREAM  at  the 
stagnation  points  and  its  independency  of freestream  mach  numbers.  ZSTREAM  is  a  finite 
element  based  streamline  code,  which  is  Mach  number  dependent  and  uniformly  valid 
everywhere  on  the  body  surface.  It  is  capable  to  define/plot  the  complete  streamline  domain  on 
body  surface,  including  the  stagnation  point,  according  to  surface  flow  solutions  given  by  a 
panel  code  for  example,  ZONA7U)  or  a  CFD  code  (For  example,  CFL3D). 


The  input  required  is  the  surface  velocities  at  each  grid/element  points.  The  output  is  the  surface 
streamlines.  ZSTREAM  functionality  is  to  provide  streamlines  input  for  Aeroheating/Heat- 
transfer  programs  such  as  SHABP/MARKV'  or  MINI  VAR  for  computations  of  the  heat-transfer 
rate  at  the  body  surface. 

Two  test  cases  are  presented  to  demonstrate  the  resulting  streamlines  and  the  developed 
aeroheating  procedure  using  ZONA7U  +  ZSTREAM  +  SHABP  for  Cp  and  heat-transfer  rate 
predictions.  These  cases  are: 

A)  3D  Validation:  Blunt  Cone  Case  and 

B)  2D  Verification:  Blunt  Wedge  Case 

Given  flow  condition  reads: 

M=10.6,  5c  =  15°  and  a  =  10°.  Solution  outputs  for  validation  are  the  pressure 
coefFicent  Cp  and  the  heat-transfer  rate,  q . 

3,1  ZSTREAM  Development 

Several  problems  were  found  in  the  QUADSTREAM  code  of  SHABP/MARKV.  First,  the  code 
only  works  for  quadrilateral  elements.  Second,  the  Newtonian  steepest  decent  method  used  in 
QUADSTREAM  calculates  streamlines  only  based  on  the  freestream  velocity  vector  and  the 
normal  vector  of  the  element.  No  local  flowfield  variables  are  involved.  We  found  the  method 
is  not  accurate,  nor  is  it  uniformly  valid  throughout  the  surfaces. 

By  contrast,  there  are  two  other  approaches  that  will  yield  more  accurate  streamline  solutions 
that  the  Newtonian  steepest  decent  method  used  in  QUADSTREAM.  One  is  to  use  the  inviscid 
surface  pressure,  and  the  other  is  to  use  the  inviscid  surface  velocity.  Both  can  yield  streamlines 
that  are  mach  number  dependent  and  uniformly  valid  throughout.  ZSTREAM  adopts  the  latter 
as  input  to  yield  higher  quality  streamline  solutions. 
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3.2  Methodology  of  Tracing  Streamlines 


A  streamline  on  a  three-dimensional  surface  is  traced  by  a  step-marching  method  using  the  flow 
velocity  on  the  tangential  plane  of  the  surface.  In  order  to  define  the  tangential  plane  at  a  point, 
the  geometry  of  the  body  surface  is  approximated  by  a  number  of  four-noded  quadrilateral 
elements  (Fig  3.1)  and  the  surface  close  to  the  nose  is  modeled  by  some  three-noded  triangular 
elements  (Fig  3.2) 


Fig  3.1  Quadrilateral  Isoparameter  Element 


Fig  3.2  Triangular  Isoparameter  Element 
Over  an  element,  coordinates  can  be  determined  using  the  following  equation 


x.(4,n)=l,NM,iX  (1) 

a=l 

where  Ne  is  the  number  of  element  nodes  (4  for  quadrilateral  elements  and  3  for  triangular 
elements),  are  the  coordinates  at  the  node  a,  ^  and  p  are  parameters  ranged  from  -1  to  1,  and 
N„  (^,  T|)  are  shape  functions  given  by 
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W=(#,'7)  =  7(1  +  I)(I-7) 

;  (2) 

W,(^,7)  =  -(l+f)(l  +  7) 

for  quadrilateral  elements,  and 

+  (3) 

for  triangular  elements. 

A  local  Cartesian  coordinate  system  x'  is  defined  on  the  element  with  x,'  along  ^  direction  and 
x'  along  the  normal  direction  (see  Fig.  3.3). 


n 


Fig.  3.3  Local  Coordinate  System 


Based  on  this  coordinate  system,  coordinates  (x,  y)  and  velocity  (u,  v)  at  the  element  nodes  on 
the  tangential  plane  can  be  calculated  in  terms  of  their  global  values  by 
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(4) 


;=i 

/=t4X 

;=i 

;=i 

M 

where  g“  are  the  global  components  of  the  velocity  and  Lij  and  L2j  are  the  direction  cosines  of 
the  local  coordinate  axes  x[  and  x'2  with  respect  to  the  global  coordinate  system. 

The  streamline  marches  from  a  beginning  point  towards  the  stagnation  point  along  the  negative 
direction  of  the  tangential  velocity  (u,  v).  Denoting  the  current  position  with  (xo,  yo)  and  the 
velocity  at  this  point  with  (uo,  vo)  and  referring  to  Fig.  3.4,  the  new  position  (x,  y)  can  be 
determined  using  the  relation; 


(^-^oK  =(t-ToK 


(6) 


Fig.  3.4  Marching  from  position  yo)  to  (x,  j) 


In  equation  (6),  the  coordinates  x  and  y  in  the  tangential  plane  can  be  calculated  using  shape 
functions  as  follows 

a=\ 


(7) 
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The  distance  between  the  current  position  {xq,  yo)  to  the  new  position  (x,  y)  is  controlled  by  step 
size  or  A7.  Using  {^0,  tjo)  to  denote  the  parameters  at  the  current  position,  their  values  at 
the  new  position  can  be  written  as 


TJ  =  7Iq+A7J 


(8) 


Thus,  for  a  given  step  size  we  can  solve  equation  (6)  for  Arj  and  vice  versa.  This  task  is 
accomplished  by  using  the  software  Mapple.  For  quadrilateral  elements,  the  results  are: 


lixi  =  [vO  (4  «x0  -  (-1  +  deta  +  etO)  trxlt  (-1  +xL0)  +  (-1  +deta  +  etO)  tx2it  (1  +  xiO)  - 
(1  +  deta  +  etO)  ★  (x3  +x4  +  x3  ★  xiO  ~  x4  frxiO))  + 
aO  fr  (-4  tryO  +  (-1  +  iieta+  etO)  t  (-1  +  xiO)  wyl  -  (-1  +  data  +  et0)  tr  (1  +  xiO)  wy2  + 

(1  +  deta-f  etO)  *  (y3  +xi0*y3  +  y4  -  xiO  *y4))]  /  [vO^ 

((-1  +  deta  +  etO)  trxl  +  x2  +  x3  -  x4  ~  (data  +  etO)  n  (x2  - x3  +x4))  +  uO  iryl  - 
aOt^  (y2  +y3  +  (deta  +  etO)  t  (yl-y2  +  y3))  +  (1  +  deta  +  etO)  *u0  wy4] 

(9) 

for  the  given  step  size  A  7  (here  denoted  by  deta),  and 


iieta=  -[vO  (-4x0+  (-1  +  etO)  xl  (-1+ dxi  +  xLO)  -(-l  +  et0)x2  (l  +  dxi  +  xLO)  + 

(1+etO)  (x4  -x4  (dxi  +xiO}  +x3  (1  +dxi  +  xiO)))  - 
^0  (-4y0  +  (-1+etO)  (-1+dxi  +xlO)  yl+  (1  +  dxi  +xiO)  (y2  -et0y2  +  y3  +et0y3)  - 
(1+  etO)  (-1  +  dxi  +  xiO)  y4)]  /  [vO  (-x2  +  x3  +  x4  +  xl  (-1+  dxi  +  xiO)  - 
(x2-x3  +x4)  (dxi  +xiO))  +u0  (-(-1  +  dxi  +  xlO)  yl+  (1  +  dxi  +  xiO)  (y2  -  y3)  + 
(-1  +  dxi  +xiO)  y4)] 

for  the  given  step  size  A(^  (here  denoted  by  dxi). 

The  results  for  a  triangular  element  are: 


dxi  =  [vOiY  (-2  itxO  +x2  -  (data  +  etO)  tr  (xl  -x3)  +  x3  +  (-xl  +  x2)  t^xiO)  + 
aO*  (2*y0  +  detatryl  +  etO  fryl +xiO  fryl-y2  -xi0*y2- 
(l+deta  +  et0)fry3)]  /  [vO  tr  (xl-x2)  +  aO  w  (-yl  +  y2)] 


(11) 


for  the  given  step  size  A  7 ,  and 

data  =  [vO  w  (-2  trxO  -dxi  w  xl-  etO  trxl  +  x2  +  dxi  trx2  +  x3  +  etO  trx3  - 
xli^xiO  +x2  ★xiO)  +  uO  w  (2  *y0  +dxi  wyl+  etO  t  yl+  xLO  *yl  - 
y2  -  dxi  IV  y2  -  xLO  *  y2  -  (1  +  etO)  w  y3)  ]  /  [vO  tv  (xl  -  x3)  +  uO  tv  (-yl  +  y3)  ] 

for  the  given  step  size  A(^  . 


(12) 
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Once  A<^  and  A;/  are  obtained,  the  new  position  can  be  determined  by  equation  (7)  and  the 
tangential  velocity  at  the  new  position  can  be  calculated  by 


a=l 

a=l 


(13) 


This  process  is  repeated  until  the  element  boundary  is  reached.  On  the  element  boundary,  ^ 
and/or  rj  equal(s)  + 1 .  Once  this  condition  is  satisfied,  the  next  nearest  element  is  considered  in  a 
similar  way,  until  the  stagnation  point,  characterized  by  u=v=0,  is  encountered.  Here  are  some 
examples; 


Fig.  3.5  Body  and  Wing  Streamlines  Produced  by  ZAERO/ZSTREAM 
(a)  Ogive-Cylinder  Body  (M=6.0,  0=2°,  t=0.022)  (b)  Simple  Wing  Body  Combination  at  qfO 


Fig.  3.6  Blunt-Nose  Cone  Streamlines  Produced  by  ZSTREAM 
(Two  different  views  with  flow  condition  at  M=10.6,  <i<=10°,  ^0=15°,  Rn=1.1‘") 
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3.3  Aeroheating  Procedure  Demonstration 

Two  cases  are  selected  for  case  study;  namely,  the  15  “-Blunt  Cone  and  the  15°-Blunt  Wedge 
with  the  given  conditions  below. 


Given  Hypersonic  Flow  Conditions  for  Blunt  Cone/Wedge 

Mach  No. 

M  =  10.6 

Wedge/Cone  angle 

II 

o 

AoA 

a=W 

Meridian  Plane 

^=180° 

Nose  Radius 

Rn=1.1” 

Freestream  Pressure 

P.  =  2.66P*^ 

Freestream  Temperature 

T„  =  89.971°^ 

Wall  Temperature 

T„  =  540”'' 

Our  objective  is  to  validate/verify  the  results  obtained  with  existing  results/test  data  as  collected 
by  Moore  (Ref  32,  See  Fig.  3.1) 

3.4  Pressure  Distributions,  Cp 

Case  A:  Blunt-Nose  Cone 

To  validate  the  computation  procedure  with  the  existing  data,  we  select  the  CFL3D  code,  for 
flow  accuracy,  to  generate  the  surface  pressure.  Fig.  3.4.2  shows  Cp  versus  surface  distance,  at 
the  meridian  plane  cp  =  180®,  up  to  18  times  of  nose  radius.  The  demarcation  points  on  the  body 
exhibit  the  following; 

•  Sonic  Point  ^  *  at  x  =  0.4  Rn 

•  Tangent  Point  at  x  =  0.74  Rn 

•  First  Minimum 

Pressure  at  x  =  2.6  Rn 

•  Sharp  Cone  Pressure 

Restored  at  x  =  12  Rn 

(Cpc=0.378) 
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Accordingly,  several  observations  can  be  made: 

i)  In  the  nose  region,  the  CFL3D  result  checks  quite  well  with  the  Newtonian-Lees 
distribution’^  (Ref  3).  One  exception  is  that  the  flow  characteristics  at  6 1  does  not  reflect 
a  kink  in  Cp. 

ii)  In  the  aft-body  downstream  an  asymptotic  matching  point  exists  (say,  X/Rn  =  12),  the  Cp 
recovers  the  value  of  that  due  to  a  sharp  cone  of  the  same  angle  Sc  (called  equivalent 
cone).  It  is  known  that  the  sharp  cone  Cp  typically  serves  as  an  asymptotic  pressure  for 
the  blunt-nose  cone  problem. 

iii)  The  relaxation  distance,  defined  as  the  distance  from  the  first  minimum  Cp  location  to  the 
asymptotic  matching  point,  is  not  known  in  advance. 

Case  B:  Blunt-Nose  Wedge 

The  present  2D  case  is  a  heuristic  study  with  the  goal  of  understanding  the  physics  of  the  blunt- 
nose  body  flow  characteristics  while  keeping  a  relatively  low  CFD  computing  level.  Further,  the 
Cp  obtained  could  be  immediately  linked  with  the  Heat-Transfer  code  (MARKV).  No  coupling 
with  the  QUADSTREAM  is  needed,  since  the  2D  body  shape  is  the  basic  streamline  itself 

Fig.  3.9  presents  Cp  versus  X/Rn  up  to  x  =  4.0  Rn,  where  Cp  converges  to  the  asymptotic 

matching  value  of  an  equivalent  wedge  of  the  same  S^,  =  15°).  The  demarcation  points  on  the 

wedge  exhibit  the  following; 

•  Sonic  Point  ^  *  at  x  =  0.47  Rn 

•  Tangent  Point  ^at  x  =  0.74  Rn 

•  First  Minimum 

Pressure  at  ^lor  x  =  0.74  Rn 

•  Sharp  Wedge  Pressure 

Restored  at  x  =  3  Rn 

(Q,.  =  0.45) 

Similarly,  observations  can  be  made; 

i)  Similar  to  the  Blunt-Cone  Case,  the  NLD-Q,  formula  should  be  applicable  to  the  wedge 
case  up  to  6. 

ii)  The  asymptotic  matching  value  in  the  aft  body  is  indeed  the  sharp  wedge  value. 

iii)  Again  the  relaxation  distance  is  better  defined  in  this  case,  as  its  starting  point  is  the 
tangent  point  (defined  by  the  body  symmetry)  6\ 


+  Also  known  as  the  improved  modified  Newtonian  approach 
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Fig.  3.9  Cpfor  1.1-in.  Nose  Radius  15°  Half-Angle  Wedge  at  M„=10.6,  a=  10°,^!?=180° 


3.5  Heat  Transfer  Rates,  q 

With  the  computed  pressures,  we  then  used  the  ZSTREAM  to  couple  them  with  the  MARKV  of 
SHABP.  Shown  in  Figs.  3.10  and  3.11  are  the  Heat  Transfer  Rates,  q,  for  a  15°-Blunt  Cone 
(Case  A)  and  for  a  15°-  Blunt  Wedge  (Case  B).  Comparisons  of  the  present  results  with  that  of 
the  existing  methods  (MINIVAR  (Ref  33),  AEROHEAT  (Ref  34))  and  Cleary’s  test  data  (Ref 
35)  are  made.  A  fairly  good  trend  is  found  between  the  present  result  and  the  test  data. 
However,  in  both  q -distributions,  large  change  in  the  heat-transfer  gradient  is  found  for  the 
present  (ZONA)  solution  at  about  one  radius  distance  from  the  nose.  The  Blunt-Wedge  (Case  B) 
Case  is  well  understood  in  that  the  large  gradient  point  corresponds  to  the  tangent  point  0  i, 
where  the  pressure  also  suffers  from  a  high  gradient  change. 

However,  for  the  Blunt-Cone  Case  (Case  A),  the  cause  of  such  a  high  gradient  in  q  is  not  clear 
which  requires  further  investigation. 
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Fig.  3.10  q  for  1.1-in.  Nose  Radius  15°  Half-Angle  Cone  at  M„=10.6,  a=  10°,  cp=  180' 


Fig.  3.11  q  for  1.1-in.  Nose  Radius  15°  Half-Angle  Wedge  at  Moo=10.6,  «=  10°,  f<=  180' 


SECTION  4 


Application  of  ZONA  Hypersonic  Aerothermodynamic  Method  to  X-34 


Summary 

•  The  developed  ZONA  Hypersonic  aerothermodynamic  method  consists  of  coupled  software 
of  ZONA  7U  +  ZSTREAM  +  SHABP.  To  demonstrate  the  developed  methodology,  we  apply 
the  software  to  X-34  for  solution  validation  with  that  of  CFD  aerothermodynamic  method 
(For  example,  CFL3D+LATCH)  in  terms  of  pressure  distributions  and  temperature 
distributions  on  the  front/lee/wind-side  surfaces  of  X-34,  at  two  selected flow  conditions. 

•  Streamlines  of  X-34  at  these  conditions  are  displayed.  Note  that  the  X-34  wing-body  adopted 
for  demonstration  is  defined  as  the  partial  configuration  (Fig  4. 1(b))  by  NASA-LaRC  (Ref  13, 
A1AA98-0880) 


a)  Full  configuration  used  for  viscous  solutions. 


b)  Partial  configuration  used  for  inviscid  solutions. 


Figure  4.1  X-34  Configurations  Taken  From  AIAA  98-0880  (Ref  13) 
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4.1  Pressure  Coefficients  of  X-34  Wing-Body  Configuration  at  M  =  6.0,  a  =  9°  and 
15.22° 

The  X-34  wing-body  configuration  consists  of  a  body  with  a  round  nose  of  7.0”  radius,  a  strake, 
and  a  swept  wing.  The  inviscid  surface  Cp  distributions  computed  by  ZONA7U  and  CFL3D 
Euler  solver  on  the  wind-side  and  lee-side  of  X-34  are  shown  in  Fig  4.2  for  the  M  =  6.0,  a  =  9° 
case  and  Fig  4.3  for  the  M  =  6.0,  a  =  15.22°  case.  In  general,  ZONA7U  results  on  the  wind-side 
correlate  well  with  the  CFL3D  results  but  those  on  the  lee-side  present  discrepancies.  This  is 
caused  by  the  flow  separation  on  the  lee  side  at  high  angle  of  attack  (a  =  9°  and  15.22°)  that 
cannot  be  handled  well  by  the  ZONA7U  attached-flow  assumption.  Nevertheless,  in  the  high 
pressure  regions  such  as  at  the  round  nose  of  the  body,  and  along  the  leading  edge  of  the  wing, 
ZONA7U  results  give  good  agreement  with  those  of  the  CFL3D.  Note  that  those  high-pressure 
regions  represent  the  high  temperature  area  on  the  configuration,  and  therefore,  are  critical  for 
the  TPS  design. 
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_ _ (a)  Front  View 

ZONA7U 


CFL3D 


^.04  0.12  028  044  06  0.76  0.92  1.08  1.24  1.4  1.56 


_ (b)  Wind-Side 

ZONA7U 


CFL3D 


■0.04  0.12  028  0.44  0.6  0.76  0.92  1.08  1,24  1.4  1.56 


■0.04  0.06  0.2  0,32  0,44  0,56  0,66  0,8  0,92  1,04  1,16  1,26  1,4  1,52  1,64 


(c)  Lee-Side 


Figure  4.2  Inviscid  Surface  Pressure  Distributions  on  tiie  X-34  at  M„  =6,  a=9*; 
(a)  Front  View,  (b)  Wind-Side,  and  (c)  Lee-Side. 
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_ (a)  Front  View 

ZONA7U 


CFL3D 


-0,04  0.12  0.28  0.44  0.6  0.76  0.92  1.08  1.24  1.4  1.56 


ZONA7U 


(b)  Wind-Side 


-0  04  012  0.28  0.44  0.6  0.76  0.92  1.06  1,24  1.4  1.56 


CFL3D 


_  (c)  Lee-Side _ _ 

Figure  4.3  Inviscid  Surface  Pressure  Distributions  on  Ae  X-34  at  =6,  a=15.22®; 
(a)  Front  View,  (b)  Wind-Side,  and  (c)  Lee-Side. 
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(a)  a=9®,  (b)  a=15.22°. 

4.2  Streamlines  and  Temperature  Distributions  on  X-34  Wind  Body:  M=6.0,  a=9° 
and  15.22® 

Fig  4.4  depicts  the  streamlines  computed  by  ZSTREAM  for  a  =  9°  and  15.22°  cases.  The 
aeroheating  analysis  is  performed  using  the  radiative  equilibrium  temperature  wall  boxmdary 
condition  (the  “hot  wall”  condition)  with  emissivity  of  0.8.  The  temperature  distributions 
computed  by  ZONA7U  +  SHABP  and  CFL3D  +  LATCH  on  the  wind-side  and  lee-side  of  X-34 
at  M  =  6.0,  are  presented  in  Fig  4.5  for  the  a  =  9°  and  the  altitude  =183  Kfl.  condition  and  in  Fig 
4.6  for  the  a  =  15.22°  and  altitude  =  112  Kft  condition.  In  both  cases,  the  temperature 
distributions  computed  by  ZONA7U  +  SHABP  on  the  wind-side  of  X-34  correlate  well  with 
CFL3D  +  LATCH  whereas  some  discrepancies  are  shown  on  the  lee-side  of  X-34. 
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ZONA7U+SHABP 


CFL3D+LATCH 


_ (a)  Front  View 

ZONA7U+SHABP 


CFL3D+LATCH 


0  100  200  300  400  500  600  700  800  900  1000  1 100  1200 


_ (b)  Wind-Side 

ZONA7U+SHABP  ^ 


CFL3D+LATCH 


(c) 


Lee-Side 


Figure  4.5  Turbulent  Surface  Temperatures  (°F)  on  the  X-34  at 
Mg,  =6,  a=9®,  Alt.=183  Kft;  (a)  Front  View,  (b)  Wind-Side,  and  (c)  Lee-Side. 


40 


ZONA7U+SHABP 


CFL3D+LATCH 


Figure  4.6  Turbulent  Surface  Temperatures  (®F)  on  the  X-34  at 
M„  =6,  a=15.22®,  Alt.=112  Kft;  (a)  Front  View,  (b)  Wind-Side,  and  (c)  Lee-Side. 
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SECTION  5 


TRAJECTORY  ANALYSIS 


Summary 

Z0NA7U  +  SHABP  is  used  to  compute  the  heat  rate  at  the  stagnation  point  of  the  X-34, 
according  to  two  assigned  trajectories  Sg,  namely  X1004601  and X1004701.  Here  ZONA7U  + 
SHABP  only  requires  the  trajectory  to  be  submitted  once,  then  it  outputs  the  pressure  Cp  and  the 
heat  rate  q  throughout  the  complete  trajectory.  For  14  time  steps  with  a  stretch  of 800  seconds, 
it  requires  10  minutes  of  computing  time.  By  contrast,  MINIVER  requires  manual  input  for  each 
point  of  interest,  i.e.  each  output  q  solution  requires  approximately  5  to  10  minutes. 


5.1  Descriptions  and  Functionality  of  Blocks  5:  Trajectory  Analysis  (Block  5) 

Description: 

-  Two  trajectory  analysis  codes  used  within  industry  (OTIS  and  POST/OTIS)  have  often 
been  adopted  by  the  Air  Force  Research  Laboratory  (AFRL).  The  trajectory  code  POST 
will  be  adopted  in  the  program  integration. 

-  POST  (Program  to  Optimize  Simulated  Trajectories,  Version  II)  is  a  generalized  point 
mass,  three  and  six-degree-of-freedom  (3-6  DOF)  trajectory  program.  The  3-DOF- 
trajectory  simulation  can  be  used  to  compute  loads,  propellant  requirements,  propulsion 
and  aerodynamic  trim  interactions  for  ascent  trajectories. 

-  POST  II  has  a  discrete  parameter  targeting  an  optimization  capability  that  can  be  used  to 
guide  the  user  to  the  optimum  trajectory  (dynamic  pressure  profile)  and  fuel  and  oxidizer 
settings  to  minimize  propellant  requirements  for  a  given  mission. 

-  POST  II  can  model  multiple  powered  and  unpowered  vehicles  including  SSTO,  TSTO, 
VTHL,  HTHL,  re-entry  problems,  as  well  as  exo-atmospheric  orbital  transfer  problems. 

Functionality: 

-  The  main  function  of  the  trajectoiy  analysis  is  to  obtain  an  optimal  trajectory  that 
minimizes  the  fuel  while  satisfying  other  constraints  such  as  Mach  number  needed  for 
specific  engine  usage,  final  velocities,  altitudes,  launch  angle,  etc.  (Ref  36). 

Inputs  include  forces,  moments  of  TRIM  aerodynamics  and  mass  from  ASTROS*.  Other 
inputs  based  on  the  requirements  from  propulsion,  mission  and  TPS  performance 
considerations. 

-  Provide  trajectory/flight  condition  history  as  inputs  to  TPS  for  heat  rate  estimate  (see  cases 
studied) 
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5.2  Phase  I  Achievements  in  Trajectory  Analysis/Case  Demonstration  Trajectory 
Analysis  Example  (Block  4) 

Here,  ZONA7U  +  SHABP  is  used  to  compute  for  the  heat  rate  at  the  stagnation  point  of  the  X-34 
according  to  two  assigned  trajectories  (X1004601  and  X1004701).  Good  correlation  is  found 
between  the  present  ZONA7U  +  SHABP  method  and  MDSTTVER  (Fig  5.1).  Clearly,  the 
advantage  of  the  ZONA7U  +  SHABP  method  over  MINTVER  is  two  fold: 

1 .  Automated  generation  of  heat  rate  solutions  once  the  trajectory  history  is  given  as  input 

2.  Formulation  of  the  aerodynamic  (Cp)  and  aerothermodynamic  {q)  is  based  on  a  rigorous 
panel  methodology  which  accounts  globally  for  the  aerod)mamics/aerothermodynamics  of 
the  complete  X-34  configuration. 


Figure  5.1  Heat  Rate  Comparison  (hot  wall)  at  Stagnation  Point 
(a)  X1004601,  (b)  XI 004701,  (c)  Trajectory  and  flight  condition  history. 
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5.3  Case  Study  (D):  Trajectory  Analysis/Case  Demonstration 

The  example  case  demonstrates  the  heat  rate  computation  methodology  in  the  presence  of  two 

given  trajectory  profiles.  The  trajectory  cases  X1004601  and  X1007401  were  selected  from  a 

previous  NASA  study  (Ref  36)  upon  the  recommendation  of  OSC  (Orbital  Sciences 

Corporation). 

The  following  is  noted  in  the  present  cases  studied: 

1.  Trajectory  X1004601  has  a  peak  Mach  number  of  7.4  and  a  maximum  AoA  at  25°, 
trajectory  X1004701  has  a  peak  Mach  number  of  8.6  and  a  maximum  AoA  at  30°. 

2.  Similar  to  the  previous  NASA  study,  we  simply  computed  the  heat-rate  at  the  stagnation 
point  of  X-34,  only  to  demonstrate  the  capability  of  the  present  methodology. 

3.  We  employ  MINIVER  and  ZONA7U  +  SHABP  for  heat-rate  computation.  Hot-wall 
condition  is  imposed,  where  the  emissivity  is  given  as  s  =  0.8. 

4.  ZONA7U  +  SHABP  only  requires  the  trajectory  inputs  to  be  submitted  once,  then  it  outputs 
the  pressure  (Cp,  not  shown)  and  the  heat-rate  (q)  solutions.  For  14  time  steps  along  a 
stretch  of  800  seconds,  it  requires  less  than  10  minutes  of  computing  time.  By  contrast, 
MINIVER  requires  manual  input  for  each  point  of  interest;  i.e.,  each  output  q  solution 
requires  approximately  five  to  ten  minutes. 
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SECTION  6 


TPS  SIZING 


Summary 

•  The  TPS  sizing  objective  is  to  develop  a  procedure  to  minimize  the  TPS  weight  while 
satisjying  the  thermal  protection  requirement  and  the  load-carrying  requirement  of  the 
combined  TA  V/TPS  structure.  An  elementary  TPS  sizing  procedure  can  be  established  by 
constructing  a  prototypical  TPS/AFRSI  model. 

•  With  this  model,  the  objective  becomes  one  to  minimize  the  weight  of  the  middle  insulation 
layer  (3).  The  inequality  constraints  are  the  operating  temperatures  of  each  layer  including 
that  of  the  skin  layer.  The  TPS  element  is  selected  on  the  windward  centerline  of  X-34.  The 
model  input  is  the  heat  rate,  which  is  provided  by  ZONA7U+SHABP  from  trajectory 
aerothermodynamic  prediction.  Based  on  a  hot-wall  consideration,  maximum  temperatures 
in  outer/interior  skin  layers  and  minimum  TPS  weight  are  resulting  outputs  obtained  by 
applying  an  iterative  procedure  using  MINIVER/EXITS.  A  formal  TPS  optimization  method 
using  Complex-Variables  Differentiation  Sensitivity  with  MINIVER/EXIT  is  proposed  in 
Section  9. 


6.1  Description  and  Functionality  of  TPS  Sizing  (Block  4) 

Description: 

-  Factors  that  influence  the  design  of  TPS  include  TPS  weight,  heat  transfer  property, 
mechanical  property,  manufacturability,  and  price. 

-  In  the  present  design  concept,  there  should  be  at  least  two  disciplines  involved  in  the  TPS 
design:  the  heat  transfer  analysis  for  TPS  size  and  the  stress  analysis  for  TPS  strength. 

-  MINIVER/EXITS  is  used  for  heat  transfer  analysis  and  ASTROS*  is  used  for  stress 
analysis. 

-  TPS  debonding  from  the  main  structure  is  of  concern.  Shear  stress  applied  on  the  TPS  and 
aeroelastic  behavior  of  the  vehicle  should  be  accurately  predicted  form  the  ASTROS* 
computation  outputs  in  order  to  prevent  the  TPS  from  a  debonding  design  under  load. 

-  TPS  is  non-sizable  with  respect  to  the  ASTROS*  main  structure  design  cycle.  It  will  be 
sized  independently;  while  from  a  structural  standpoint,  it  will  be  considered  jointly  with 
the  entire  structure. 

Functionality: 

-  Objective  is  to  minimize  the  TPS  weight  while  satisfying  the  thermal  protection 
requirement  and  the  load-carrying  requirement  of  the  combined  TAV/TPS  structure. 

-  Provide  feedback  to  ASTROS*  for  the  load  carrying  design:  TPS  mass,  stiffness,  material 
degradation,  temperature  mapping  to  the  main  structure. 

-  Shock/Shear  loads  input  from  ASTROS*  for  TPS  heat/stress  analyses. 

-  Heat  rate  input  from  trajectory  analysis  (Block  5) 
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6.2  Phase  I  Achievements  in  TPS  Sizing  Case  Demonstration:  TPS  Sizing  Example 
(Block  4) 

1.  In  order  to  facilitate  the  case  studied,  ZONA  has  integrated  the  EXITS  module  of  MINIVER 
with  ZONA7U+SHABP.  In  so  doing,  this  automated  procedure  allows  the  heat  rate 
solutions  of  the  latter  to  be  directly  input  into  EXITS  for  heat  transfer  analysis. 

2.  The  following  TPS  sizing  example  (case  studied)  demonstrates  that  a  TPS  sizing  procedure 
has  been  developed  and  tested  out  successfully  on  a  typical  case  (TPS  element)  previously 
adopted  by  NASA  (Ref  37).  This  case  yielded  a  minimum  weight  insulated  layer,  among  the 
three  such  layers  selected,  resulting  in  realistic  thermal  solution  in  the  structure  layer  (skin). 
A  flow  chart  is  attached  showing  such  a  optimization  concept  utilizing  EXITS  for 
TPS/AFRSI  sizing  (see  Fig  6.1) 

3.  The  case  studied  suggests  that  present  procedure  can  be  generalized  into  the  following 

it  poses  no  restriction  on  the  number  of  selected  insulated  layers  to  be  considered 
the  thickness  of  the  other  layers  of  the  TPS  can  be  adjusted  at  will 
it  should  also  include  shear-stress  consideration 

it  can  be  developed  into  an  automated  optimization  scheme  for  TPS  minimum  weight 
objective  with  required  thermal/mechanical  constraints. 


Figure  6.1  Trial-and-Error  Procedure  to  Obtain  an  Optimized  AFRSI. 
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6.3  Case  Study  (C):  TPS  Sizing  Case  Demonstration 


Example  Case:  TPS  of  X-34 

Case  example  of  a  TPS  unit  is  selected  from  Ref  37  (pg  7  NASA  TM  2000-210289).  The 
AFRSI  concept  modeled  in  the  example  is  composed  of  an  outer  fabric  with  C-9  coating,  6  Ib/ft 
Q-fiber  felt  insulation,  and  an  inner  fabric  layer,  and  it  is  attached  to  the  structure  with  RTV 
adhesive  (Fig  6.2). 

This  TPS  unit  is  placed  at  the  windward  bottom  surface  of  the  X-34  centerline  at  50  in.  from  the 
nose  (see  Fig  6.3,  Pt.  A).  The  X-34  is  subject  to  a  trajectory  with  heat  rate  limiting  at  point  A. 

Layer  1  -  Coating  (0.01  in.  HRSI  Coating) 

Layer  2  -  Outer  Fabric  (0.015  in.  Outer  Fabric  AB3 12) 

Layer  (3)  Insulation 

a.  Q-Felt  Insulation  (standard) 

b.  Q-Felt  3.5PCF 

c.  6LB  Dynaflex 

(Insulation  layer  size  is  to  be  optimized) 

Layer  4  -  Inner  Fabric  (0.009  in.  Inner  Fabric  AB312) 

Layer  5  -  Adhesive  (0.008  in.  RTV  Adhesive) 

Layer  6  -  Structure  (0.011  in.  Aluminum) 

.  Tauter  and  T interior  are  the  temperatures  at  the  outer  edge  and  (1)  to  (5)  interior  layers  of  the  TPS. 
Tskin  Is  the  temperature  at  the  nodes  within  the  skin  layer  6.  _ _ 

Figure  6.2  Description  of  the  AFRSI  on  the  Structure. 


t 

X  (inches) 

i 


Centerline  (bottom  view  and  side  view,  L=50  in.) 
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The  TPS/AFRSI  Model 

According  to  the  above,  a  TPS/AFRSI  model  is  constructed  in  Fig  6.2.  The  modeled  AFRSI 
consists  of  six  layers  with  layer  three  being  the  insulation  layer  subjected  to  adjustment  and  layer 
six  being  the  structure  (or  the  skin)  layer,  within  which  a  temperature  is  set  (Tupper-bound  —  300  F) 
as  a  constraint. 

The  Iteration  Process 

In  the  present  case,  we  fix  all  other  layers;  the  objective  then  becomes  the  minimum  weight  of 
layer  three,  the  insulation  layer.  The  inequality  constraints  are  the  operating  temperatures  of 
each  layer  including  the  skin  layer  (Tskin  ~  300°  F).  Each  layer  has  a  Tupper-bound  specific  to  the 
material  composition. 


Discussion  of  Results 

In  Table  6. 1 ,  the  thickness  of  insulation  materials  used  as  layer  3  is  adjusted  to  have  maximum  of 
Tskin  within  0.1%  of  300°  F  (the  constraint  of  the  skin  material.  Aluminum)  throughout  the  entire 
trajectory.  The  initial  temperature  is  set  as  100°  F  and  the  inner  edge  is  assumed  insulated. 


Table  6.1  Thickness  and  Weight  Comparison  Between 


Layer  3 
material 

Thickness 

Normalized 
weight,  TPS 

Normalized 
weight,  layer  3 

Max 

Touter 

Max 

T  interior 

Max 

Tskin 

Q-Felt  insulation 

0.456  in 

1.000 

1.000 

708.7°  F 

696.4°  F 

300.3°  F 

Q-Felt3.5PCF 

0.638  in 

0.694 

0.408 

713.6°  F 

702.0°  F 

300.2°  F 

6LB  Dynaflex 

0.560  in 

1.118 

1.228 

696.9°  F 

r  mT\  on 

681.6°  F^ 

300.2°  F 

(Tcutcr  and  Ti„,erior  are  the  temperatures  at  the  outer  edge  and  interior  oj  ll':>.  l  skin  ts  me  temperature  at  ine 
nodes  within  the  skin.  “Maximum  temperature’’  is  determined  by  scanning  all  temperatures  obtained 

throughout  the  trajectory  history.)  •  •  •  u 

Computation  results  indicate  that  Touter  and  Tjnterior  reach  maximum  at  t  —  340  sec,  which  is  right 
after  the  peak  heating  condition  of  the  XI 004601  trajectory  at  t  =  330  sec.  (ZONA7U  +  SHABP 
result  has  been  verified  at  t  =  340  sec.) 
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Because  of  AFRSI,  the  Tskin  reaches  its  maximum  at  t  =  680  sec.  The  temperature  drop  from  the 
outer  edge  of  AFRSI  to  skin  is  about  400  °F. 

Table  6.1  shows  that  an  iteration  process  can  be  used  to  find  a  AFRSI  design  (using  Q-Felt 
35PCF  as  insulation  layer)  which  outperforms  current  AFRSI  design  (using  Q-Felt)  with  69%  of 
original  weight  from  the  heat  transfer  analysis.  This  is  a  trial-and-error  procedure  and  is  a  time 
consuming  process.  In  Phase  II,  we  will  develop  an  optimization  driver  for  MINIVER/EXITS 
that  can  automatically  achieve  a  minimum  weight  design  of  the  TPS  system. 
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SECTION  7 


TRIM  ANALYSIS 


Summary 

To  demonstrate  one  of  the  multifunctional  capabilities  of  ASTROS*,  we  used  its  Trim  Module  to 
perform  trim  analysis  of  a  modeled  X-34,  at  a  selected flight  condition  (M—6. 0,  AoA  =9°, 
h=183Kft)  With  X34/FEM  model  provided  by  OSC/NASA-Langley,  ASTROSVTRIM  results  in  a 
trim  condition  requiring  trailing  edge  flap  angle  3=^2. 05°,  load factor  N=0.97g,  and  total  weight 
of  vehicle  W=16, 000  lbs.  The  aerodynamic  loads  at  the  trim  condition  is  then  mapped  to  the 
FEM  grid  by  means  ofASTR0S*/3D  Spline  Module,  thus  allowing  a  subsequent  analysis  of  the 
X-34  structures.  Accordingly,  the  resulting  stress  distribution  of  the  X-34  at  time  is  obtained. 


7.1  Descriptions  and  Functionality  of  Blocks  3:  ASTROS*  Structural  Optimization 
Description: 

-  As  shown  in  Fig  7.1  below,  ASTROS*  consists  of  at  least  8  comprehensive  modules; 
these  include  the  Structural  FEM  module,  the  ZAERO  Aerodynamic  module,  the 
Aeroelastic  Stability  module,  the  ASE  module,  the  Optimization  module,  the  Trim 
module,  the  Smart  structures  module,  and  the  Sensitivity  module. 

-  The  outcome  of  the  proposed  ASTROS*  optimization  is  a  minimum- weight  structural 
design. 

-  It  uses  the  Aerodynamic  Influence  Coefficient  (AIC)  matrices  from  ZONA7U/AIC  to 
obtain  steady/unsteady  airloads  and  it  performs  Trim  analysis  to  yield  flight  loads. 

-  It  performs  aeroelastic  analysis  to  yield  Flutter/Divergence  constraints. 

It  will  receive  the  shock  impingement  location  from  a  high  fidelity  Unified  hypersonic 
aerodynamics. 


50 


Smart  Structures  Module 


Modeling  of  P2T  and  SMA 
activations  and  computes 
the  induced  aerodynamic 
control  forces. 


ZAERO  Aerodynamic 
Module 


Unified  steady/unsteady 
aerodynamics  for  subsonic, 
transonic,  supersonic,  and 
hypersonic  flows. 


Sensitivity  Analysis 


Provides  sensitivity  of 
stresses,  stability,  and 
performance  with 
respect  to  structural 
design  variables. 


Structural  Finite  Element  Module 
NASTRAN^ompatible  FEM  analysis 


Smart  structures  module 
^  y 


ASE  Module 


State-space  aeroservo elastic  analysis  including  SISO/MIMO 
control  system  for  stability  and  gust  loads  analysis. 


Trim  Module 

Static  aeroelastic 
analysis  for  flight 
loads  and  optimum 
trim  solutions 


Aeroelastic  Stability 
_ Module _ 

Provides  true  damping 
flutter  solutions. 


Optimization  Module 

An  optimizer  driving  all 
modules  to  achieve  the 
optimum  design. 


Figure  7.1  ASTROS*  Engineering  Modules. 


Functionality: 

-  The  flight  conditions  for  defining  the  aeroelastic  constraints  are  provided  by  the  results  of 
the  trajectory  analysis. 

-  For  TAV/TPS  design,  the  effects  of  the  TPS  mass  and  the  material  property  de^adations 
and  the  temperature  on  the  load-carrying  structure  due  to  aero-heating  must  be  included  in 
the  FEM  model. 

-  The  objective  fimction  of  the  ASTROS*  structural  optimization  is  the  minimum  weight, 
which  is  a  feedback  information  to  the  trajectory  analysis. 

-  With  the  aerod5mamics  input  from  Block  1,  the  flight  loads  including  the  shear  loads  and 
shock  loads  computed  by  the  trim  module  of  ASTROS*  is  another  feedback  information 
to  the  stress  analysis  of  the  TPS  design. 

-  The  total  mass  of  this  design  is  the  feedback  information  to  the  trajectory  analysis  for  the 
next  design  iteration.  Other  important  ASTROS*  results  are  the  loads  and  structural 
deformation  computed  by  the  ASTROS*/Trim  module  for  the  TPS  stress  analysis. 

Several  trim  solutions  computed  by  the  trajectory  analysis  will  be  selected  as  the  flight 
conditions  for  defining  the  aeroelastic  constraints.  These  constraints  include  the  flutter 
and  divergence  at  various  Mach  number  and  the  strength  constraints  of  each  designed 
structural  element  at  various  angles  of  attack  and  Mach  numbers. 


51 


7.2  Phase  I  Achievements  in  X-34  TRIM  Analysis  using  ASTROS*/FEM:  Block  3 

The  objective  of  this  study  is  to  demonstrate  the  trim  analysis  capability  of  ASTROS*  for  the 
flexible  X-34  in  hypersonic  flow.  The  outcomes  of  the  trim  analysis  are  the  control  surface 
deflection  angles,  load  factors,  etc.  as  well  as  the  stress  distribution  in  the  structures.  Fig  7.2 
depicts  an  ASTROS*  finite  element  model  of  the  X-34  that  was  converted  from  a 
MSC/NASTRAN  X-34  model  originally  provided  by  the  Orbital  Sciences  Corporation.  Using 
the  aerodynamic  influence  coefficient  (AIC)  matrix  computed  by  ZONA7U  in  conjunction  with 
the  stiffness  matrix  of  the  FEM  model,  the  ASTROS*  trim  analysis  shows  that  in  order  to  trim 
the  X-34  at  M  =  6.0,  a  =  9°  and  altitude  =  183  Kft,  the  required  trailing  edge  flap  angle  is  2.05° 
degrees  and  a  load  factor  of  0.97-g  for  a  total  weight  of  16,000  lbs.  At  this  condition,  the 
aerodynamic  loads  computed  at  the  ZONA7U  panels  are  then  mapped  to  the  FEM  grid  using  the 
3D  spline  module;  allowing  a  subsequent  stress  analysis  of  the  structures.  Such  a  stress 
distribution  is  shown  in  Fig  7.3. 
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SECTION  8 


CONCLUDING  REMARKS  AND  RECOMMENDATION  FOR  PHASE  II 

WORK 

1.  Central  Methodologies  (Blocks  1-5)  required  for  aerothermodynamic  optimization  were 
individually  developed  and  validated. 

2.  Hypersonic  Aerodynamics  and  Aerothermodynamic  methodology  for  Blocksl/2  are  aimed  at 
replacing  the  high-level  method  CFL3D+LATCH  and  the  low-level  method  MINIVER  by  a 
mid-level  method ZONA7U+  SHABP  which  has  been  developed  in  Phase  I. 

3.  For  ZONA7U+SHABP  to  generate  one  set  of  X-34  aerodynamic/heat  rates  typically  requires 
10  minutes  on  a  550  MHZ  PC,  whereas  for  CFL3D+LATCH  it  requires  30  hours. 

4.  ZSTREAM  was  developed  in  Phase  I  to  replace  the  QUADSTREAM  streamline  generator  of 
SHABP  in  that  the  latter,  derived  from  the  Newtonian  flow  consideration,  has  a  stagnation- 
point  singularity  in  its  streamline  solution  and  it  does  not  depend  on  freestream  Mach 
number. 

5.  Validation/verification  of  ZONA’s  aerothermodynamic  method  ZONA7U+SHABP  suggests 
that  further  improvement  is  needed  in  the  following:  ZONA7U  requires  higher-fidelity 
upgrade  in  order  to  cope  with  the  high  AoA  and  the  lee-side  aerodynamics  of  SHABP  needs 
to  be  replaced  by  tht  AEROHEAT methodology  (Ref  34)  in  order  to  further  improve  the  local 
heat  rate  estimates. 

6.  The  TPS  weight  sizing  example  shows  that  the  designed  TPS  weight  can  be  further  reduced 
if  an  automated  optimized  scheme  can  be  developed.  A  database  of  TPS  material  in  terms  of 
their  thermal  and  mechanical  properties  must  be  fully  established  in  order  to  enhance  the 
capability  of  the  optimized  scheme. 

7.  The  trim  solution  of  the  X-34  in  terms  of  the  flight  loads,  input  to  the  structural  FEM  within 
ASTROS*,  will  yield  shear  loads  and  shock  loads  which  will  result  in  strength  constraint  in 
the  ASTROS*  optimization  procedure. 

8.  Given  trajectory  inputs,  ZONA7U+SHABP  aeroheating  solution  at  the  nose  of  X-34  was 
verified  with  previous  solutions  obtained  by  NASA  (Ref  36).  Total  optimization  loop 
including  ASTROS*  will  be  tested  next  using  an  X-34  example  as  a  demonstration  case. 
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SECTION  9 


PHASE  II  PLAN 


Summary 

In  Phase  /,  the  ZONA  unified  hypersonic/supersonic  aerodynamic  method  ZONA  7U  in  the 
aerothermoelastic  software  development  for  TPS/TA  V  design/analysis,  was  proven  a  successful 
tool  through  feasibility  cases  studied  in  the  CKEM  body,  blunt  cones  and  X-34  wing-body. 
Further  upgrade  of  ZONA  7U  to  a  high-fidelity  panel  method  ZONAIR  is  proposed  in  order  to 
enhance  the  universality  of  the  aerodynamic  capability  with:  1)  unified  AIC  for 
hypersonic/supersonic/subsonic  Mach  numbers;  2)  exact  paneling  to  model  OML  surface;  2) 
high  AoA  flow;  4)  corrected  location/strength  of  impinged  shocks;  5)  two-body  interference 
aerodynamics.  Other  proposed  improvements  of  the  aerothermoelastic  software  include;  i) 
incorporation  ofAERJOHEA  T  to  replace  SHABP  with  ZONAIR  and  ZSTREAM;  ii)  optimized  TPS 
sizing  procedure  using  Complex  Variable  Differentiation  Sensitivity  with  MINIVER/EXIT;  Hi) 
applying  temperature  mapping  onto  structural  FEM.  The  structural  optimizer  ASTROS*^  is  used 
throughout  for  sizing  TA  V/TPS  while  maintaining  the  TPS  design  as  a  part  of  the  load-carry 
structure.  The  TPS  design  adopts  POST/OTIS  for  trajectory  optimization  in  the  outmost  design 
cycle,  thereby  leading  to  a  final  minimized fuel  and  weight  objective.  Meanwhile,  ZONAIR  and 
ASTORS*  FEM  will  be  driven  by  the  automated  parametric  giid-panel  mesh  generation  of  AML, 
developed  by  Technosoft.  By  means  of  AML,  an  integrated  software  framework  with  a  feature- 
based  environment  can  be  established  to  support  an  underlying  object-oriented  architecture  that 
will  harness  and  link  all  software  modules  in  configuration  geometry,  aerothermoelastic 
analysis,  TPS  design,  etc.  RMLS  and  HyperX  are  proposed  as  the  test-bed  vehicles  for  software 
demonstration  in  Phase  II. 


9.1  Phase  II  Technical  Objectives 
Main  Goal: 

The  ZONA  team  proposes  to  continue  the  phase  I  development  in  a  unified 
hypersonic/supersonic  aerothermoelastic  methodology,  to  enhance  its  geometry  high-fidelity  and 
to  integrate  it  with  a  TAV/TPS  structural  design/optimization  procedure  for  TAV/TSP  weight 
minimization  throughout  TAV’s  re-entry/maneuver  flight  phases. 

Further,  the  above  development  will  be  integrated  in  a  feature-based  design  environment  with 
parametric  control  of  models  and  data  exchange  using  the  Adaptive  Modeling  Language  (AML). 

Specific  Objectives: 

1.  Development  of  an  expedient  conceptual/preliminary  design  tool  that  allows  varying  of 
parametric  geometry  for  rapid  assessment  of  design  concepts. 
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2.  The  TPS  sizing  is  to  meet  the  thermal  protection  requirements  of  the  main  structures  while 
being  included  with  the  main  structures  and  treated  together  as  a  part  of  the  load-carrying 
structure. 

3.  The  developed  software  can  be  utilized  by  a  design  engineer  as  an  outer  loop  design 
procedure  for  trajectory  optimization,  which  in  turn  will  lead  to  a  minimum  fuel  and  weight 
TAV  design. 


9.2  Phase  II  Work  Plan 

The  ZONA  team  proposes  a  two-year  effort  with  eight  tasks  for  Phase  II.  A  block  diagram 
showing  a  work  plan  for  performing  these  eight  tasks  is  depicted  in  Fig  9.1. 


Figure  9.1  Block  Diagram  of  the  Phase  II  Work  Plan. 


•  AML 

In  Phase  II,  a  framework  with  a  feature-based  design  environment  will  be  developed  supporting 
an  underlying  object-oriented  architecture  that  will  link  various  aspects  of  vehicle  configuration, 
geometry,  analysis,  and  sizing  using  the  Adaptive  Modeling  Language  (AML)  developed  by 
TechnoSoft,  Inc. 

The  design  environment  will  support  the  parametric  design  of  the  vehicle  outer  mold  line  and 
primary  internal  structure  geometry.  Fully  automated  parametric  grid-panel-based  mesh 
generation  for  aerodynamic  analysis  and  structured  meshes  for  structure  analysis  will  be 
integrated. 

•  TPS  Sizing 

The  environment  will  provide  links  to  various  analysis  applications  through  a  common 
computational  model.  These  links  will  support  the  automated  exchange  of  information  among 
the  analyses,  geometry  and  grids.  Links  to  different  analysis  tools  will  be  provided  including; 
aerodynamic,  aeroheating,  trajectory,  thermal  protection  system  (TPS)  sizing,  and  structural 
sizing  (all  included  in  the  Block  diagrams  of  Fig  1.1  and  Fig  9.1).  A  closure  model  supporting 
the  direct  link  among  the  analysis  tools  and  iteration  control  will  be  incorporated  for  vehicle 
sizing. 
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Based  on  an  initial  design  configuration,  aerodynamic  loads  are  computed  and  fed  into  a 
trajectory  analysis.  The  results  of  the  trajectory  analysis  are  used  as  input  for  TPS  analysis  and 
sizing  algorithms.  The  TPS  configuration  will  be  computed  at  selected  points  on  the  vehicle 
outer  mold  line  (OML). 

While  considered  jointly  with  the  TPS  as  a  load-carrying  structure,  the  main  structure  will  be 
optimized  subjected  to  strength  requirements  based  on  a  set  of  critical  load  conditions.  After  the 
TPS  and  main  structure  are  determined,  an  updated  vehicle  weight  is  created.  Based  on  this  new 
weight  model,  a  new  trajectory  analysis  is  initiated.  This  automated  loop  iterates  until  the  TPS 
design  converges. 

•  ZONAIR,  ZSTREAM,  andAEROHEAT 

The  aerothermodynamic  analysis  is  performed  using  the  ZONAIR  code,  a  high-fidelity,  high- 
order,  unified  subsonic/supersonic  panel  method.  In  Phase  II,  we  will  enhance  ZONAIR  to 
include  hypersonic  aerodynamics  by  incorporating  the  hypersonic  methodology  of  ZONA7U 
into  ZONAIR.  This  enhancement  is  discussed  in  section  4.1.  The  enhanced  ZONAIR  will  be 
integrated  with  ZSTREAM  to  provide  inviscid  surface  pressure  distribution  and  streamlines  to 
the  AEROHEAT  code,  a  more  accurate  aeroheating  tool  than  SHABP.  The  background  of 
AEROHEAT  is  discussed  in  section  4.2.  The  integrated  ZONAIR  +  ZSTREAM  +  AEROHEAT 
code  will  be  used  to  generate  an  aerodynamic  and  temperature  database  for  trajectory  analysis 
performed  either  by  POST  or  OTIS.  The  resulting  time  history  of  the  aeroheating  data  computed 
by  POST/OTIS  will  be  the  input  of  the  MINIVER/EXITS  code  for  TPS  sizing  at  various  vehicle 
locations.  In  Phase  II,  we  will  develop  an  optimization  driver  of  MINIVER/EXITS  to  automate 
the  TPS  sizing  procedure  using  an  innovative  complex-variable  differentiation  sensitivity.  The 
formulation  of  the  complex-variable  differentiation  scheme  will  be  presented  in  section  4.3 

After  the  design  is  completed,  the  TPS  stiffness  and  mass  effects  will  be  included  in  the  FEM 
model  using  the  approach  discussed  in  section  4.4.  This  FEM  model  is  generated  using  the 
automated  parametric  grid-panel-based  mesh  generator  of  AML  discussed  in  section  4.5. 
Meanwhile,  the  temperature  distribution  resulting  from  the  TPS  design  will  be  mapped  onto  the 
surface  elements  of  the  FEM  model.  This  mapping  procedure  is  fully  automated  using  a 
procedure  presented  in  section  4.6. 

•  ASTROS^ 

ASTROS*  will  be  used  to  perform  the  structural  optimization  of  the  main  structures.  However, 
prior  to  this  optimization  computation,  an  ASTROS*/trim  analysis  will  first  be  performed  at 
various  conditions  along  the  vehicle  trajectory  to  obtain  the  critical  load  conditions.  Strength 
constraints  at  these  critical  load  conditions  along  with  the  other  aeroelastic  instability  constraints 
will  be  imposed  to  achieve  an  optimum  structural  design.  The  results  of  the  structural  design  and 
TPS  sizing  will  allow  for  the  creation  of  an  updated  weight  model.  The  analysis  closure  loop  is 
then  re-initiated  with  this  new  weight  model  and  the  process  iterates  until  convergence  is 
reached. 
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•  Test  Beds 

To  validate  the  proposed  design  environment  and  demonstrate  its  rapid  design  capability,  we 
propose  to  apply  the  developed  software  on  two  configurations,  namely,  the  RMLS  and  the 
Hyper-X  as  test  beds.  Section  9.8  presents  twelve  subtasks  outlining  a  step-by-step  design 
procedure  to  conduct  this  task.  The  outcome  of  this  task  will  be  a  preliminary  RMLS  design. 


9.3  High  Fidelity,  High-Order  Panel  Method  for  Hypersonic/Supersonic 
Aerodynamics 

As  shown  in  the  Phase  I  study  of  the  15°  blunt  cone  and  X-34  cases,  ZONA7U  loses  its  accuracy 
at  high  angles-of-attack  due  to  the  attached-flow  formulation.  With  continuous  R&D  in  panel 
method  development  since  1995,  ZONA  has  developed  a  unique  software  product  called 
ZONAIR,  a  high  fidelity  unified  subsonic/supersonic  panel  code.  ZONAIR  has  the  following 
capabilities. 

•  Unstructured  Grids  for  Modeling  Complex  Air  Vehicles  such  as  Conventional  Aircraft 
Blended  Wing-Body  and  TAV Configuration 

ZONAIR’ s  paneling  scheme  is  based  on  the  so-called  “unstructured”  grids  as  opposed  to  the 
“structured”  grids  adopted  by  other  high-fidelity  panel  codes  such  as  PANAJR.  Fig  9.2  presents 
the  comparison  between  ZONAIR’ s  unstructured  paneling  scheme  and  PANAIR’ s  paneling 
scheme  where  the  advantages  of  adopting  the  unstructured  grids  are  described. 


ZONAIR 

Unstructured  Grids 

CQUAIM 

PANAIR 

Structured  Grids 

Wiflg/betfy  tAtty-wt  «aka  _  / 

Body  (MM  ^  y'  y 

v— Piwmo  bas* 

•  Similar  to  structural  FEM  (MSC/NASTRAN),  the 
entire  configuration  is  defined  by  “grids”. 
CTRIA3’s  and  CQUAD4’s  define  the  connectivity 
between  the  grids. 

•  Only  the  starting  lines  of  the  wake  need  to  be 
defined  (via  GEAR  elements).  There  are  no  input 
requirements  for  the  surface  wake. 

•  PATRAN,  FEMAP,  etc.,  can  be  employed  directly 
for  pre-  and  post-processing.  ZONAIR  outputs 
PLOAD4  bulk  data  cards  for  pressures,  local  Mach 
numbers,  and  velocities  on  each  aerodynamic  panel. 

•  The  entire  configuration  is  first  divided  into  several 
“networks”.  Each  network  is  further  divided  by  m  x 
n  set  of  grids. 

•  The  location  of  the  wake  surfaces  must  be  explicitly 
defined. 

•  No  commercially  available  software  can  be  used 
directly  for  pre-  and  post  processing. 

Figure  9.2  •  Comparison  of  ZONA! 

R  and  PANAIR  Paneling  Schemes. 
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Figure  9.3  Regular  and  Random  Paneling  of  a  Sphere  at  M=0.0  and  a=0.0  deg. 


Another  advantage  in  using  unstructured  grids  is  that  it  allows  arbitrary  grid  point  selection  for  a 
given  configuration.  In  order  to  demonstrate  this  feature,  a  sphere  is  modeled  by  using  regularly 
spaced/shaped  panels  (called  Regular  Panels)  and  randomly  spaced/shaped  panels  (called 
Random  Panels)  whose  pressure  distribution  results  are  shown  in  Figs  9.3(a)  and  9.3(b), 
respectively.  Clearly,  this  arbitrary  grid  point  selection  capability  of  the  unstructured  grids  can 
greatly  reduce  the  user  burden  in  the  grid  generation  process. 

•  Vortex  Roll-Up  for  High  Angles-of-Attack  Aerodynamics 


FREE  PARAMETERS:  12.  (3.  I4,  (g,  Yc  2^ 

_ FIXED  PARAMETERS!  62.  63,  64.  % _ 

9.4  Free/Fed  Vortex  Sheet  Kinematics  for  Vortex  RoU-up. 


Fig  9.4  shows  the  essential  elements  in  ZONAIR  for  vortex  roll-up  modeling;  the  vortex  sheet 
emerging  from  the  wing  leading  edge  and  tip  (free  vortex  sheet)  and  the  rolled-up  core  or  spiral 
region  (fed  vortex  sheet)  fed  by  the  leading-edge  and  tip-vortex  sheets.  The  following  boundary 
conditions  are  imposed  on  these  elements: 


-  The  configuration  surface  must  be  impermeable. 

The  free  sheet  and  wake  cannot  support  a  pressure  difference  and  must  be  impermeable  as 
well. 

-  The  fed  sheet  is  an  extension  of  the  free  sheet  and  feeds  vorticity  to  the  vortex  core  (modeled 
as  a  simple  line  vortex).  The  boundary  condition  governing  the  fed  sheet  size  and  core 
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orientation  is  that  the  total  force  induced  on  the  fed  sheet  and  core  by  the  rest  of  the 
configuration  be  parallel  to  the  core. 

-  Kutta  conditions  are  imposed  along  the  appropriate  leading,  side,  and  trailing  edges  of  the 
wing  in  the  presence  of  free  sheets  emanating  from  these  edges. 

Using  the  above  vortex  roll-up  model,  the  induced  drag,  lift,  and  pitch  moment  versus  angles-of- 
attack  of  a  generic  advanced  fighter  (GAF)  at  M  =  1.8  is  shown  in  Fig  4.1.4.  It  can  be  seen  that 
ZONAIR’s  prediction  correlates  very  well  with  the  wind  tunnel  data  up  to  a  =  20°. 

•  Unified  Hypersonic/Supersonic  ZONAIR 

Currently  ZONAIR  is  based  on  the  potential  flow  equation,  which  is  valid  only  in  the  supersonic 
Mach  number  range.  In  Phase  H,  we  will  incorporate  the  hypersonic  methodology  of  ZONA7U 
into  ZONAIR,  that  is: 

•  Equivalent  Mach  number  transformation  to  circumvent  the  superinclined  panel  problem. 

•  Local  pulsating  body  analogy  to  include  the  flow  rotationality  effects. 

Additionally,  ZONAIR  will  be  further  improved  on  its  hypersonic  capability  in  two  aspects  by 
using  a  Perturbed-Euler  formulation  (Ref  21). 

•  Correct  the  impinged  shock  strength  and  location  (e.g.,  nose  shock  of  the  X-34  impinging 
on  the  wing). 

•  Accoimt  for  the  cross-flow  near-normal  shock  effect  on  the  lee-side  due  to  high  AOA 
flow  (in  addition  to  the  vortex  roll-up  effect). 

Shown  in  Fig  9.5  are  two  such  PEF  shock  solutions  versus  that  of  CFL3D.  Without  ofrier 
secondary  shock  interaction,  the  PEF  could  correct  the  nose  shock  on  the  wing  to  provide  more 
accurate  shock  load. 


Figure  9.5  Solutions  of  Perturbed  Euler  Formulation  (PEF)  vs  CFL3D, 

(a)  Concave  Body  at  M=50,  y=1.4, 6=38°,  and  (b)  Convex  Body  at  M=15, 7=1.4, 6=35°. 


Once  the  above  hypersonic  methodology  is  incorporated  into  ZONAIR,  the  three  test  cases 
shown  in  section  2.0,  namely  the  CKEM  body,  the  15°  blunt  cone,  and  the  X-34  wing-body 
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configuration  will  be  employed  to  validate  ZONAIR  at  high  angles-of-attack  and  for  the 
correction  of  impinged  shock  location. 

The  enhanced  ZONAIR,  once  developed,  can  serve  as  a  universal  aerodynamic  tool  for 
aerodynamic  force/moment  generation,  streamline  generation  for  aeroheating  analysis,  AIC 
generation  for  aeroelastic  analysis,  and  trim  solutions  for  critical  loads.  Its  high-fidelity  paneling 
scheme  is  compatible  with  FEM  mesh  generation,  whereas  its  computing  time  is  two-orders  less 
than  that  of  CFD  codes.  Table  9. 1  shows  the  comparison  of  functionality  between  ZONAIR  and 
other  aerodynamic  codes  where  by  comparison  ZONAIR’ s  universality  and  superiority  is  clearly 
seen. 


A  ZONAIR  (Wave  Drag) 

O  Wind  Tunnel  Data  (Total  Drag) 


ANGLE  OF  ATTACK 


Figure  9.6  Force  and  Moment  Coefficients  of  GAF  vs  Angle-of- Attack  at  M  —  1.8 


Table  9.1  Comparison  between  ZONAIR  and  Other  Aerodynamic  Codes. 


Code 

Method 

Compu¬ 

tational 

Efficiency 

Streamline 

Solution 

For 

Aeroheating 

Hypersonic/ 

Supersonic/ 

Subsonic 
Mach  No, 

AIC  for 
Structural 
FEM 

Geometry 

High 

FideUty 

High 

AOA 

2  Body 
Aero 
Inter¬ 
ference 

CFL3D 

Euler/N-S 

30  hrs/ 
X-34 

“■  Yes 

All 

No 

Yes 

Yes 

Yes 

PAN  AIR 

Potential 

20  mia^ 
X-34 

No 

Supersonic/ 

Subsonic 

No 

Yes 

No 

Yes 

ZONAIR 

Potential  + 
PEF 

20  min/^ 
X-34 

Yes 

All 

Yes 

Linear- 
Order  Panel 

Yes 

Yes 

ZONA7U 

Potential  + 
PEF 

10  min/  I 
X-34 

Yes 

All 

Yes 

Constant- 
Order  Panel 

No 

Yes 

APAS 

Potential  +  ! 
Empirical 

<  10  min 

Newtonian 

S.L. 

Empirical  for 
hypersonics 

No 

Low-Order 

Panel 

No 

Yes 
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MINIVER 

Analytical/ 

Empirical 

«  10  min 

No 

No  subsonics 

No 

No 

No 

No 

DATCOM 

Analytical/ 

Empirical 

«  10  min 

No 

All 

No 

No 

Yes 

No 

AP98 

Analytical/ 

Empirical 

«  10  min 

No 

All 

No 

No 

Yes 

No 

9.4  Adopting  the  AEROHEAT  Code  for  Accurate  Aerothermodynamic  Analysis 

SHABP’s  empirical  equations  for  aeroheating  analysis  are  based  on  two  simple  hypersonic 
similarity  solutions;  one  for  stagnation  point  and  one  for  flat  plates.  As  shown  in  our  Phase  I 
study  on  the  CKEM  body  and  15°  blunt  cone,  these  similarity  solutions  are  not  accurate 
particularly  near  the  stagnation  point. 

In  Phase  II,  we  will  replace  SHABP  by  a  more  accurate  aeroheating  code  called  AEROHEAT. 
AEROHEAT  solves  the  convective-heating  equations  using  an  axisymmetric  analogy  that  allows 
any  axisymmetric  boundary  layer  method  to  be  applied  along  an  inviscid  surface  streamline. 
These  convective-heating  equations  consist  of  a  set  of  approximate  convective-heating  equations 
developed  by  Zobyl  (Ref  38)  that  provides  accurate  surface  heating  rates  with  a  minimal  amount 
of  computational  effort.  Laminar  and  turbulent  heating  rates  are  calculated  by  relating  the 
surface  skin  friction  to  the  momentum  thickness  Reynolds  number.  Note  that  the  3D  effects  of 
the  AEROHEAT  methodology  is  included  through  the  streamline  metric  coefficients,  which  can 
be  accurately  provided  by  the  ZSTREAM  code.  The  inviscid  aerodynamic  solutions  required  by 
AEROHEAT  will  be  computed  using  ZONAIR.  The  integrated  ZONAIR  +  ZSTREAM  + 
AEROHEAT  solutions  will  be  validated  with  the  CFL3D  +  LATCH  results  and  test  data  for 
CKEM  body,  15°  blunt  cone  and  the  X-34  wing-body  configuration. 

9.5  Optimization  Procedure  for  TPS  Sizing 

In  Phase  I  we  conducted  a  TPS  sizing  study  on  the  X-34  using  the  EXITS  module  of  MINIVER 
(as  described  in  section  6).  Based  on  the  temperature  constraint  of  300°  F  on  the  skin  layer,  we 
have  shown  that  a  minimum  weight  TPS  can  be  achieved  by  a  trial-and-error  procedure. 
However,  performing  such  a  trial-and-error  procedure  for  minimum  weight  TPS  design  over  the 
whole  vehicle  is  a  very  time  consuming  process.  In  Phase  II,  we  will  develop  an  automated  TPS 
design  procedure  using  the  feasible  direction  method  as  an  optimization  driver  of 
MINIVER/EXITS  to  achieve  the  optimum  design. 

Fig  9.7  depicts  a  typical  TPS  sizing  problem 
that  consists  of  n  layers  with  thickness  hi, 
hj—hn,  respectively.  For  a  given  heat  flux  q 
applied  on  the  outer  boundary,  the  objective  is 
to  minimize  the  total  weight  of  the  TPS  system 
while  keeping  the  temperature  at  each  layer 
(7))  below  their  respective  maximum 
operational  temperature,  lou  The 

corresponding  optimization  formulation  is 
shown  as  follows; 


j  Layer  1 

Layer  2 

X 

k 

/ 

h„ 

/ 

Layer  n 

_ To _ 

Figure  9.7  Typical  TPS  Sizing  Problem. 
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Minimize:  JV  =  where  /?,  is  the  density  of  the  layer 

/=1 

Subjected  to:  T,  <  Tot  i  =  1,  2... n. 

Design  variables:  /?,  >  0  i  =  1,  2...  n. 

This  optimization  problem  can  be  effectively  solved  by  the  feasible  direction  method  if  the 
sensitivity  of  the  constraint  function  can  be  evaluated  accurately  and  efficiently.  The 
MINIVER/EXITS  module  can  be  used  to  compute  the  temperature  distribution  at  each  layer,  but 
it  can  not  provide  the  sensitivity  of  the  temperature  distribution  with  respect  to  the  design 
variables  hi. 

To  obtain  this  sensitivity  information,  we  will  apply  the  complex-variable  differentiation 
approach  on  the  MINIVER/EXITS  module.  The  complex-variable  differentiation  technique, 
first  proposed  by  Lyness  and  Moler  (Ref  39)  and  extended  to  BEM  solution  recently  by  ZONA 
(Ref  40),  is  a  powerful  mathematical  technique  that  can  provide  the  “numerically  exact” 
derivatives  of  a  complicated  function.  In  this  complex  variable  approach,  the  variable  /?  of  a  real 
function  T(h)  is  replaced  by  a  complex  one,  h  +  iAh.  For  a  small  Ah,  T(h+Ah)  can  be  expanded 
into  a  Taylor’s  series  as  follows: 

T(h  +  iAh)^T(h)  +  iAh^  +  ... 

oh 


The  first  derivative  of  the  above  function  can  be  expressed  as: 

dT  _Im(T(h  +  iAh))  I  , 
dh  Ah 

It  can  be  seen  that  the  derivative  using  the  complex  variable  approach  only  requires  function 
evaluation  whereas  the  conventional  finite  differencing  technique  involves  the  differencing  of 
two  functions.  Therefore,  the  complex  variable  approach  does  not  suffer  from  the  cancellation 
errors,  hence,  becomes  step-size  independent  in  the  small  asymptotic  liiuit.  In  fact,  the  step  size 
can  be  as  small  as  machine  zero  (10‘^^)  and  still  retain  accuracy  (Ref  40). 

To  incorporate  the  complex  variable  technique  into  the  MINIVER/EXITS  module  for  sensitivity 
analysis  is  straightforward  simply  by  declaring  all  variables  in  the  MINIVER/EXITS  module  as 
complex  variables.  The  imaginary  part  of  the  thickness  input  of  MINIVER/EXITS  represents  a 
small  incremental  thickness  whereas  the  sensitivity  is  the  imaginary  part  of  the  temperature 
output  divided  by  the  incremental  thickness. 

9.6  Inclusion  of  TPS  Mass  and  Stiffness  Effects  in  FEM  for  Optimization  of  Main 
Structures 

As  discussed  earlier,  the  TPS  design  is  primarily  based  on  the  thermo-protection  constraints,  not 
the  structural  constraints.  Detailed  structural  design/analysis  of  the  TPS  system  can  be 
performed  only  if  the  load  paths  from  the  main  structures  to  the  TPS  are  known.  Therefore,  such 
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a  structural  design/analysis  of  TPS  should  be  performed  in  the  detailed  design  stage.  On  the 
other  hand,  because  the  TPS  is  placed  on  the  vehicle  outer  mold  line  that  gives  significant  area 
moment  of  inertia,  excluding  the  TPS  stiffhess  from  the  FEM  model  can  lead  to  a  over-designed 
main  load-carry  structure. 

To  overcome  the  above  issue,  in  Phase  II  we  will  develop  a  procedure  that  can  automatically 
convert  the  overall  stiffiiess  and  mass  of  TPS  into  equivalent  composite  laminate  properties. 
Those  laminates  will  be  included  as  the  top  layer  of  the  vehicle  composite  skin  or  in  terms  of 
ASTROS*  input  bulk  data  cards,  the  top  layer  in  the  PCOMP  bilk  data  cards  of  the  skin 
elements.  While  the  thickness  of  each  layer  in  the  skin  elements  can  be  defined  as  design 
variables  of  the  ASTROS*  structural  optimization,  these  top  layers  remain  non-sizable.  In  doing 
so,  the  TPS  structures  are  not  part  of  the  ASTROS*  optimization  but  their  load-carry  capability 
can  influence  the  size  of  the  main  load-carry  structures. 

After  the  ASTROS*  optimization  is  completed,  the  load  paths  from  the  TPS  including  the 
interlaminar  shear  stresses  (from  ASTROS*  stress  recovery  calculations)  and  the  shock  loads 
(from  the  ZONAIR  pressure  distributions  according  to  corrected  shock  location)  will  be  output 
for  the  subsequent  detailed  design  of  the  TPS  structures. 

9.7  Automated  Parametric  Mesh  Generation  for  ASTROS*  and  ZONAIR 

Automated  parametric  mesh  generation  for  MSC/NASTRAN  models  is  an  existing  capability  in 
the  Supersonic  Hypersonic  Vehicle  Design  (SHVD)  system  developed  by  LMCO/TechnoSoft. 
Because  of  the  similarity  between  the  MSC/NASTRAN  and  ASTROS*  bulk  data  cards,  a  mesh 
generator  for  ASTROS*  FEM  models  can  be  developed  with  minor  modification  to  the  SHVD 
system.  The  mesh  generator  of  ASTROS*  can  also  directly  be  employed  for  the  ZONAIR 
aerodynamic  model  generation  because  of  the  unstructured  paneling  scheme  of  ZONAIR. 

In  addition  to  the  mesh  generation,  the  ASTROS*  optimization  requires  the  selection  of  design 
variables  and  the  definition  of  constraint  functions.  The  development  of  an  automated  procedure 
to  establish  such  an  optimization  input  will  be  a  major  effort  of  this  task.  In  order  to  monitor  the 
progress  of  the  ASTROS*  optimization  computation,  we  will  develop  a  real-time  graphical 
capability  to  display  the  design  variables  on  the  FEM  model  along  with  the  active  constraints  at 
each  optimization  iteration.  If  an  optimum  solution  can  not  be  achieved,  the  graphical  capability 
will  help  the  user  to  quickly  identify  the  source  of  the  problem  and  consequently  modify  the 
optimization  problem  statement  until  an  optimum  solution  is  obtained. 

9.8  Temperature  and  Aeroloads  Mapping  from  Aerodynamic  to  Structural  Grids 

For  the  present  aerothermoelastic  analysis,  two  types  of  data  mapping  between  the  aerodynamic 
grid  (the  aerodynamic  panels)  and  the  structural  finite  element  (FEM)  grid  are  required.  The  first 
type  is  the  mapping  of  the  aerodynamic  forces  from  the  aerodynamic  grid  to  the  structural  grid  as 
well  as  the  displacement  from  the  FEM  grid  back  to  the  aerodynamic  grid.  This  type  of  data 
mapping  procedure  has  been  fiilly  developed  in  the  ASTROS*  code  that  contains  four  spline 
methods,  namely  the  infinite  plate  spline  method  (Ref  41),  the  thin  plate  spline  method  (Ref  42), 
the  beam  spline  method  and  the  rigid  body  attachment  method  that  jointly  generate  a  spline 
matrix  for  displacement  mapping  from  FEM  grid  to  aerodynamic  grid.  Based  on  the  principle  of 
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virtual  work,  the  transposed  spline  matrix  can  be  employed  for  the  force  mapping  from  the 
aerodynamic  grid  to  the  FEM  grid. 

The  second  type  mapping  requirement  is  the  temperature  mapping  from  the  aerodynamic  grid  to 
the  FEM  surface  grid.  In  Phase  II,  we  will  develop  a  finite-element-based  mapping  procedure  as 
shown  in  Fig  9.8.  This  procedure  assumes  a  bilinear  temperature  distribution  over  the 
aerodynamic  quadrilateral  panel  and  temperature  are  defined  at  the  corners  of  the  panel.  Fig  9,8 
indicates  the  equations  used  in  the  temperature  projection  process.  These  comer  temperatures  to 
the  individual  aerodynamic  panels  are  then  mapped  to  the  surface  of  the  FEM  model  that  is 
comprised  of  plate  or  membrane  elements. 

The  output  of  this  temperature  mapping  procedure  will  be  a  set  of  TEMP  and  TEMPD 
NASTRAN  or  ASTROS*  bulk  data  cards  to  define  the  temperature  at  FEM  grid  points  for 
determination  of  thermal  loading,  temperature-dependent  material  properties,  or  stress  recovery. 


Figure  9.8  Temperature  Mapping  from  Aerodynamic  to  Structure  Grids. 

9.9  Test  Beds  of  the  Proposed  Design  Environment 

The  proposed  design  environment,  once  developed,  will  be  validated  and  tested  on  two  TAV 
configurations,  a  Reusable  Military  Launch  System  (RMLS)  shown  in  Fig  9.9  and  the  Hyper-X 
shown  in  Fig  9.10.  The  Hyper-X  is  an  existing  TAV  design  whose  aerodynamic  and  structural 
models  are  shown  in  Fig  9.11.  Since  the  models  exist,  no  mesh  generation  effort  will  be 
required.  In  addition  to  the  X-34,  the  selection  of  Hyper-X  is  primarily  for  the  validation  of 
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ZONAIR+ZSTREAM+AEROHEAT  with  CFL3D+LATCH  results.  One  such  a  CFDL3D 
pressure  distribution  on  the  Hyper-X  is  shown  in  Fig  9. 12. 

The  RMLS  system  is  a  two-stage  reusable  launch  vehicle  concept  whose  exact  outer  mold  line  is 
presently  undetermined  and  is  subject  to  the  design  of  several  geometric  controlling  parameters 
such  as  fuselage  fineness  ratio,  nose  fineness  ratio,  wing  span,  wing  chord,  etc.  In  particulm, 
ZONAER  is  the  only  successful  software  that  has  a  proven  capability  in  2-body  separation  by  its 
store  separation  option  in  validation  with  wind  tunnel  and  drop  test  data  at  subsonic  and 
supersonic  speeds.  Here,  this  project  provides  an  ideal  test  bed  of  the  proposed  design 
environment  to  demonstrate  its  rapid  design  capability  for  an  optimum  and  converged  solution. 


Figure  9.9  A  Reusable  Military  Launch  System  (RMLS). 


Figure  9.10  Hyper-X  Configuration. 
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Figure  9.12  CFDL3D  pressure  distribution  on  the  Hyper-X. 


First,  based  on  a  set  of  initially  guessed  values  of  the  geometric  controlling  parameters,  a 
ZONAIR  model  and  an  ASTROS*  FEM  model  are  generated  using  the  automated  parametric 
mesh  generator.  Consequently,  an  initial  vehicle  weight  is  calculated  using  the  weight  model. 
Meanwhile,  an  aerodynamic  and  aeroheating  database  is  generated  by  performing  the 
ZONAIR+ZSTREAM+AEROHEAT  computations  at  a  set  of  Mach  numbers  and  angles  of 
attack.  This  provides  the  input  database  of  the  POST/OTIS  code  to  acquire  a  trajectory  for  the 
initial  design.  The  output  of  POST/OTIS  will  be  the  heat  rate  time  histories  that  are  used  for  the 
TPS  sizing  over  the  entire  vehicle.  The  stiffhess  and  mass  effects  of  the  TPS  system  will  be 
included  in  the  ASTROS*  main  structure  optimization.  An  updated  weight  is  then  obtained  from 
the  ASTROS*  optimization  result  and  a  new  trajectory  analysis  is  initiated.  This  process  is 
repeated  until  a  converged  solution  is  achieved. 

Next,  a  second  set  of  geometry  controlling  parameters  will  be  defined,  leading  to  a  second 
converged  solution  by  following  the  same  procedure  described  above.  The  comparison  of 
weight  and  performance  between  the  first  and  the  second  converged  solutions  will  be  conducted. 
This  vvdll  provide  the  sensitivity  of  the  RMLS  performance  with  respect  to  the  geometry 
controlling  parameters  and  provide  guidelines  regarding  changes  of  the  geometry  controlling 
parameters  to  yield  the  final  optimum  solution. 

Finally,  all  intermediate  solutions  during  the  iteration  phase  of  the  two  converged  solutions  will 
be  graphically  documented.  Based  on  experience  gained  during  the  RMLS  design,  operational 
guidelines  of  the  proposed  design  environment  will  be  established  and  final  adjustments  of  the 
software  system  will  be  made. 
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9.10  Phase  II  Statement  of  Work 


The  tasks  to  be  performed  in  Phase  II  are  defined  below.  A  schedule  of  the  tasks  is  shown  in 

section  9.11,  which  also  “shows  the  timing  for  related  program  deliverables,  meetings,  and 

presentations. 

Task  1:  Enhancement  of  ZONAIR  for  Hypersonic  Aerodynamics 

•  Incorporate  the  hypersonic  methodology  of  ZONA7U  into  ZONAIR. 

•  Validate  the  enhanced  ZONAIR  with  the  CFL3D  results  of  the  X-34  and  Hyper-X 
configurations. 

Task  2:  Integration  of  AEROHEA  T  into  ZONAIR 

•  Modify  the  AEROHEAT  code  to  accept  streamlines  and  inviscid  pressures  computed  by 
ZSTREAM  and  ZONAIR,  respectively. 

•  Validate  the  ZONAIR+ZSTREAM+AEROHEAT  results  with  CFL3D+LATCH  aeroheating 
results  of  the  CKEM  body,  15°  blunt  cone,  X-34,  and  Hyper-X. 

Task  3:  Development  of  Optimization  for  TPS  Sizing 

•  Incorporate  the  complex-variable  differentiation  sensitivity  into  the  MINIVER/EXITS 
module. 

•  Develop  an  optimization  driver  using  the  feasible  direction  method  for  TPS  optimization. 

Task  4:  Inclusion  of  TPS  Mass  and  Stiffness  Effects  for  Optimization  of  Main  Structures 

•  Develop  an  automated  procedure  to  convert  the  TPS  overall  stiffness  and  mass  into 
equivalent  composite  laminate  properties. 

•  Extract  the  inter-laminar  shear  stresses  and  shock  loads  from  ASTROS*  for  detailed  design 
of  TPS  structures. 

Task  5:  Development  of  the  Automated  Parametric  Mesh  Generator  for  ASTROS*  and 
ZONAIR 

•  Modify  the  SHVD  system  for  ASTROS*  and  ZONAIR  mesh  generations. 

•  Establish  an  automated  procedure  for  defining  the  design  variables  and  constraint  functions. 

•  Develop  a  real-time  graphical  capability  to  display  the  progress  of  the  ASTROS* 
optimization  computations. 

Task  6:  Development  of  Temperature  Mapping  Capability  from  Aerodynamic  to  Structural 
Grids 

•  Establish  a  projection  process  to  map  the  temperature  distribution  from  the  aerodynamic  to 
structural  grids. 

•  Output  the  mapping  results  in  terms  of  NASTRAN  or  ASTROS*  TEMP  and  TEMPD  bulk 
data  cards. 

Tasks  7:  Demonstration  of  the  Rapid  Design  Capability  of  the  Proposed  Design  Environment 
Using  the  Reusable  Military  Launch  System  (RMLS)  as  a  Test  Bed 

•  Generate  an  ASTROS*  and  ZONAIR  model  based  on  a  set  of  initially  guessed  values  of  the 
geometric  controlling  parameters. 
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•  Perform  a  ZONAIR+ZSTREAM+AEROHEAT  analysis  to  establish  the  aerodynamic  and 
aeroheating  database. 

•  Obtain  a  trajectory  of  the  initial  design  using  POST/OTIS. 

•  Conduct  a  TPS  sizing  over  the  entire  vehicle  based  on  the  heat  rate  time  histories  computed 
by  POST/OTIS. 

•  Perform  an  ASTROS*  optimization  computation  for  an  optimum  structural  design. 

•  Establish  an  updated  weight  of  the  vehicle  and  initiate  a  new  trajectory  analysis. 

•  Repeat  the  above  process  until  a  converged  solution  is  achieved. 

•  Obtain  a  converged  solution  by  defining  a  second  set  of  geometric  controlling  parameters. 

•  Compare  the  weight  and  performance  of  the  two  solutions  for  the  sensitivity  of  the  RMLS 
performance  with  respect  to  the  geometric  controlling  parameters. 

•  Graphically  document  all  intermediate  solutions  during  the  iteration  phase  of  the  above  two 
converged  solutions. 

•  Establish  operational  guidelines  of  the  proposed  design  environment. 

•  Conduct  a  final  adjustment  of  the  software  system  based  on  the  experience  gained  during  the 
RMLS  design. 

Task  8:  Documentation  of  the  Proposed  Design  Environment 

•  Four  manuals  will  be  generated.  The  User’s  Manual  will  include  the  input  instructions, 
modeling  guidelines,  etc.  The  complete  theoretical  formulation  will  be  presented  in  the 
Theoretical  Manual.  Sample  test  cases  for  demonstrations  of  problem  set-up  will  be 
documented  in  the  Applications  Manual.  The  Programmer’s  Manual  will  present  the  total 
program  architecture,  the  contents  of  database  data  entities,  and  the  functionality  of  each 
submodule. 
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9.11  Planned  Program  Schedule 


Tasks 


Yr  1  Quarter 
1  I  2  I  3  I  4 


Yr  2  Quarter 


12  3 


Performed 

by 


1.  Enhancement  of  ZONAIRfor  Hypersonic  Aerodynamics 

•  Incorporate  the  hypersonic  methodology  into  ZONAIR 

•  Validate  the  enhanced  ZONAIR  with  the  CFL3D  results . 

2.  Integration  of  AEROHEAT  into  ZONAIR 

•  Modify  the  AEROHEAT  code  to  accept  streamlines  and 
inviscid  pressures  computed  by  ZSTREAM  and  ZONAIR 

•  Validate  the  ZONAIR+ZSTREAM+ AEROHEAT  results . 

with  CFL3D+LATCH  aeroheating  results 

3.  Development  of  Optimization  for  TPS  Sizing 

•  Incorporate  the  complex-variable  differentiation  sensitivity 

•  Develop  an  optimization  driver  for  TPS  optimization . 

4.  Inclusion  of  TPS  Mass  and  Stiffness  Effects  for 
Optimization  of  Main  Structures 

•  Develop  an  automated  procedure  to  convert  TPS  stiffness/ 
mass  into  equivalent  composite  laminate  properties 

•  Extract  the  inter-laminar  shear  stresses  and  shock  loads . 

from  ASTROS*  for  detailed  design  of  TPS  structures 

5.  Development  of  the  Automated  Parametric  Mesh 
Generator  for  ASTROS*  and  ZONAIR 

•  Modify  the  SHVD  system  for  ASTROS*/ZONAIR  meshes 

•  Establish  an  automated  procedure  for  defining  the  design 
variables  and  constraint  lections. 

•  Develop  a  real-time  graphical  capability  to  display  the . 

progress  of  the  ASTROS*  optimization  computations. 

6.  Development  of  Temperature  Mapping  Capability  from 
Aerodynamic  to  Structural  Grids 

•  Establish  a  projection  process  to  map  the  temperature 
distribution  from  the  aerodynamic  to  structural  grids 

•  Output  the  mapping  results  in  terms  of  NASTRAN  or . 

ASTROS*  TEMP  and  TEMPD  bulk  data  cards 

7.  Demonstration  of  the  Rapid  Design  Capability  Using  the 
RMLS  as  a  Test  Bed 

•  Generate  ASTROS*/ZONAIR  model  of  the  initial  design . 

•  Perform  a  ZONAIR+ZSTREAM+AEROHEAT  analysis  . 

•  Obtain  a  trajectory  of  the  initial  design  using  POST/OTIS  ■ 

•  Conduct  a  TPS  sizing  over  the  entire  vehicle . 

•  Perform  an  ASTROS*  optimization  computation . 

•  Establish  an  updated  weight  of  the  vehicle . 

•  Repeat  the  above  process  for  a  converged  solution . 

•  .  Obtain  a  converged  solution  by  defining  a  second  set  of . 

geometric  controlling  parameters 

•  Compare  the  weight  and  performance  of  the  two  solutions  •  ■ 
for  the  sensitivity  of  the  RMLS  performance 

•  Graphically  document  all  intermediate  solutions . 

•  Establish  operational  guidelines . 

•  Conduct  a  final  adjustment  of  the  software  system . 

8.  Documentation  of  the  Proposed  Design  Environment 
Kick-Off  Meeting 

Final  Oral  Presentation  _ 


Z  =  ZONA  Technology;  T  =  TechnoSoft;  D  =  Dr.  DeJamette 
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9.12  Phase  II  Endorsement 


For  assessment  of  ZONA’s  unified  hypersonic/supersonic  unsteady  aerodynamics,  ZONA7U 
capability  the  reader  is  recommended  to  contact  the  following  users: 

-  Warren  Smith  /  Lockheed  Martin/LMMFC 
(972)  603-7726,  smith@vs.lmco.com 

-  Rudy  Yurkovich  /  Boeing  Company 

(314)  233-2563,  rudolph.yurkovich@mw.boeing.com 

-  Bart  Fowler  /  NASA-Marshall 

(256)  544-2691,  bart.fowler@msfc.nasa.gov 

Endorsement  letter  from  Dr.  Zarda  of  LMCO/LMMFC  and  from  Dr.  Henri  Fuhrmann  of  Orbital 
(OSC)  for  the  proposed  effort  are  attached  below. 

Lockheed  Martin  Endorsement  Letter _ 

Some  Excerpts... 

Lockheed  Martin  Missiles  and  Fire  Control 
supports  the  proposal  effort  ...  AFOl-252, 

Aerothermoelastic  Optimization  Methods  for 
Reusable  Launch  Vehicles. 

We  regard  built-in,  Multi-Disciplinary 
Optimization  (MD)  technology  as  a  next  generation 
capability  for  improving  air  vehicle  design  and 
manufacturing  technology. 

The  potential  impact  of  the  proposed  work 
represents  a  paradigm  shift  in  the  design  and 
corresponding  manufacturing  processes,  enabling 
greatly  improved  design  optimization. 

This  success  will  give  rise  to  a  Phase  II SBIR  effort, 
which,  if  awarded,  will  receive  our  complete 
support.  Our  organization  within  Lockheed  Martin 
anxiously  anticipates  participation  within  the 
development  and  implementation  of  the  resulting 
SA-FSV  integrated  framework  architecture. 


itCMaa 

Mhiiki  aad  fire  Cosieret  -  Orlwdo 

Mr.  Adel  Chfifftaly 
TeehnoSoft,  Inc. 

4434  Carver  Woods  Drive 
Cincinoati,  Ohio  45242 

To  whon  it  nay  eo&oem: 

Lockbeed  Manin  Misitles  and  Fire  Control  {idly  toppora:  tho  prcpotal  effort  of  TechnoSoft.  Inc., 
submioed  to  the  Uciiad  States  Airforce,  b  response  to  dw  SrsaJl  BusJoaxs  Jasovatioa  Research 
USAF  Topic  AFOl'252,  Aerothemoel^c  O^nuutjo&  Methods  for  Reusable  Lunch  Vehicles. 

As  rauaccr  of  the  En^eering  hfothods  Croup  at  Lockhw/I  Marria  Missiles  aad  Fire  Control 
(LMMAFC),  I  lead  a  team  of  englneerioB  design  experu.  These  estpern  develop  sod  udliie 
sophisticated  cagtoeehBg  tools  to  solve  complex  pf^lemi  related  to  desigslag  tad  {nasufseturing 
reissilei,  gimbsls,  teekeU,  sod  space  systems  in  geoeral.  Cortcfli  Spacx  Access  and  Fonn  Strike 
Vehicles  (SA-FSV)  Analysis  indkeuc  that  cobstantial  acrobeaiiag  is  prcscot.  AercthennalelAstjc 
nnalysU  nodThcnnal  Proiectioo  System  (TPS)  siring  of  this  pb^meooo  collrimmlvely 
contribute  to  vehiclo  eooGguniUoa  and  itructunl  layout.  Ad^tiooally,  coosidexatioa  of  the 
irarions  stage  (light  conditions  iloog  the  vehicle’s  tmjectosy  shoold  be  Indoded.  AsiduneateJ  to 
this  uDpcDvcmcai  in  design  methods  on  the  methods  that  enable  the  captere  of  design  intencrioos 
with  optimisations  within  the  imcgnicd  Adaptive  ModeTiDg  Language  (AML)  Cramewodc.  Wt 
regurt  built-in.  Multi-DiscipUniry  Optiroixaiiem  (MDO)  uxhaolo0  as  a  next  generetioo 
ca^Hity  for  improving  aJr  vehicle  design  and  taanu&ctnring  lecboology.  The  potential  impact 
of  the  proposed  work  represents  a  paradigm  shift  in  the  design  and  conespoadiog  maoufaeturiag 
processes,  enabling  greatly  Improved  design  opUiiuzatioa. 

Over  the  last  four  years,  LMM&FC  has  collahorated  with  TechnoSoft,  lac.  on  levocnl 
developmental  efforts  ecatered  on  TeehnoSoft’s  AML  Web-based  Design  Faviioimjeni  (WDE). 
Our  telationsbip  has  resulted  in  revolwiofuuy  ionovaiiotis  wjtWn  our  donuin  of  missile  de&gn. 
Our  past  experience,  coupled  with  their  sapbistJeJUad  leehnology  and  innovative  organiiatioas  of 
experu.  bas  us  convinced  that  the  Phase  I  SBIR  proposal  effort  submitted  by  TeetooSoft,  Inc.  win 
be  ensured  success.  This  success  will  give  rise  to  a  Phase  H  SBIR  effort,  whl^.  If  awarded,  will 
rceeivo  our  contplete  support  Our  orgaaizntioo  within  Lockheed  Martin  nwdously  notldpates 
pairidpatioo  within  the  devcloprocnt  and  implemenution  of  tho  rcrolting  SA-FSV  Integrated 
fcamewotk  nrchitecuue- 


Dr.  Richard  Znida 

Maonger,  Etigmeeiing  Methods  Croup 
Lockheed  Martin  Missilei  nod  Hro  Control.  Orlando 
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Orbital  Sciences  Corporation  Endorsement  Letter 


Orl^al 

Some  Excerpts... 

Novrtnl>ef  7. 2001 

Df.  D.  J  iu 

ZOMA  'IVcftnotogy 

7'!  30  H  5teKiin  Privi> 

ScotWnIc,  AZ  W251-3540 

Orbital  Sciences  Corp.  is  pleased  to  extend 
our  support  for  ZONA''  AFRL/SBIR  phase  II 
proposal . . . 

SUBJECT;  Crt>Hnl  Support  of  ZONA  Flwse  H  T’ropojal  wtitled  'OutCRinti'rf  Hypwvunlc 

A«rt»lhefjnwlas<i<:  Mothodology  for  TAV/TPS  jtnifturol  iX‘slgii  «ncl  Opluwj/ntlon." 

Denr  Or.  T.hi: 

(.'Irt’ilal  Scimros  Coiyomficr*  is  please  to  ex  lend  our  uuppart  for  ZONA 's  AFRL/SniR  phase  II  propiwftl 
cntUletJ,  "Integrated  Hypersonic  Aer(*thcrmoel"sttc  MeOnxloJogy  for  TAV/TPS  ft  fucturnt  Design  and 
OptbniraUon."  Your  work  5n  Uris  is  serving  toWdge  lire  gap  bdweiin  nnnlytical  niethixls  and  tnore 

Wnri^intensivi?  coiuputallonal  nppmnelws  for  hyfrewmic  vetiivle  desigiv  SpecificaHy,  your  product  will 
have  reUicancy  Inrywrd  the  AHlt  need  wrd  w^jll  address  a  irquiivnicnl  in  the  launch  community 
esprfialJy  within  NASA's  Spaix;  Launch  Jniliallve. 

Wc  are  pScased  that  you  liavc  been  ab!t'  to  a^etlata  gejwmtcd  lor  the  X-34  prtigram  in  the  development  of 
your  mclhixloiogy.  ATA  successfully  applied  ZONA  developed  iitethods  to  the  X.34  progrnnr during 
acroclastir  analysis  i*/  the  airframe  and  your  supjJort  of  a  troublesome  epee<}  brake  "Inirr."  preiliv.tion  was 
appreciated.  T  houglr  the  X-34  has  yet  to  bet.‘<wnc  the  na  iion«il  hypersonic  teslbcd  that  il  was  irvlcnded  for, 
much  useful  {.*xp8rim«mtn)  and  coiupnfntional  data  exists  for  just  such  efforts, 

Tlirtok  you  for  iiJh»wing  lus  the  opjjorhmity  to  view  and  comment  o»  your  work.  I  hop*  the  TPS  links  wc 
pffwidcri  have  bjion  helpful.  PleAsr  consider  u.s  as  a  resource  for  Kcuaable  L.aunch  Vehicle  and 
a«?rothermal  inforniation. 

Your  work  in  this  area  is  serving  to  bridge  the 
gap  between  analytical  methods  and  more 
labor-intensive  computational  approaches  for 
hypersonic  vehicle  design. 

Specifically,  your  product  will  have  relevancy 
beyond  the  AFRL  need  and  will  address  a 
requirement  in  the  launch  community 
especially  within  NASA's  Space  Launch 
Initiative. 

SuKvicly  /O 

Henri  D,  hullmipnn 

Flight  Dynanwfx  Manager 

Advanc^  Programs  (troup 

ATA  successfully  applied  ZONA  developed 
methods  to  the  X-34  program  during 
aeroelastic  analysis  of  the  airframe  and  your 
support  of  a  troublesome  speed  brake  "buzz" 
prediction  was  appreciated. 

For  further  information  of  ZONA  Technology  software,  please  visit  www.zonatech.com. 
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A.l  Introduction 


The  design  of  Trans  Atmospheric  vehicles  is  a  complex  multidisciplinary  problem.  The  process 
is  collaborative  and  often  requires  interaction  between  numerous  engineers  to  refine  a  design 
from  the  conceptual  to  the  preliminary  stage.  It  is  a  highly  iterative  and  interactive  process,  often 
involving  many  proprietary  engineering  applications.  This  process  involves  geometric  modeling, 
aerodynamics,  propulsion,  thermal  analysis,  trajectory  analysis,  structural  analysis,  subsystem 
layout,  and  internal  and  external  interface  specifications.  Historically,  these  analyses  and 
simulations  have  not  been  integrated  and  the  overall  design  process  has  been  inefficient  and 
poorly  defined.  TechnoSoft  is  working  with  ZONA  Technology  on  developing  an  environment 
that  will  mitigate  these  drawbacks  by  possessing  the  following  attributes; 

•  Integrated  geometry,  meshing,  aerodynamics,  aerothermal,  thermal  protection  system 
design,  weight  modules,  trajectoiy  simulation,  closure,  and  design  exploration  and 
optimization 

•  Multiple  (arbitrary  body)  aerodynamic  and  aerothermal  codes  that  range  from  subsonic  to 
hypersonic  regimes 

•  Rocket  and  air-breathing  propulsion 

•  Common  (arbitrary)  geometry  models  for  all  analysis  tools 

•  A  proven  integration  approach  for  geometry  and  analysis  tools 

•  Built-in  design  exploration  methods  that  allow  the  user  to  conduct  HAV  trade  studies  and 
multidisciplinary  optimizations 

•  An  object-oriented  architecture  that  allows  for  the  efficient  capture  of  knowledge  and 
processes 

•  An  open  architecture  that  allows  for  efficient  integration  of  legacy  application  codes 

•  Platform  independence 

•  A  web-enabled  distributed  architecture 

•  A  design  methodology  architecture  that  can  simultaneously  support  conceptual, 
preliminary,  and  detailed  designs  and  simulations  within  the  same  environment 

•  An  integrated,  feature-based  design  environment  with  a  parametric  link  to  Pro/E  that  will 
allow  geometiy  and  packaging  parameters  to  be  integrated  efficiently  with  performance 
and  cost  studies 
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\‘elucle  Design  Process 

An  architecture  and  functional  specification  were  developed  for  a  collaborative  environment 
based  on  an  underlying  object-oriented,  web-enabled,  multidisciplinary,  distributed  framework 
supporting  the  design  and  analysis  of  Trans  Atmospheric  Vehicles  (TAV).  The  developed 
architecture  and  the  prototype  system  are  tied  to  ZONA’s  hypersonic  program. 

The  developed  framework  seamlessly  integrates  several  domain-specific  analyses  and  general- 
purpose  modules  into  one  coherent  environment  that  is  linked  to  ZONA’s  analysis  codes.  The 
common  computational  model  will  allow  users  to  utilize  the  full  capabilities  of  ZONA’s  analysis 
codes  from  a  user-friendly,  web-enabled,  collaborative  environment. 

The  environment  provides  a  customizable  graphical  user  interface  supporting  a  feature-based 
design  environment  integrating  ZONA’s  analysis  codes,  and  geometry  enabling  the  rapid 
prototyping  of  TAVs.  It  supports  a  unified  geometric  part  model,  providing  various  levels  of 
modeling  fidelity  to  capture  conceptual  and  preliminary  design  processes.  The  environment 
links  multiple  users  in  a  collaborative  process,  automating  and  managing  data  transfer  and 
interactions  among  users,  designs,  analyses,  and  tools.  It  provides  multidisciplinary  optimization 
capabilities  to  enhance  vehicle  analysis,  reducing  engineering  time  and  cost  while  expanding  the 
design  space  explored.  A  common  computational  model  seamlessly  integrates  geometry  and 
ZONA’s  analysis  codes  to  support  closure  of  the  process  through  iterative  control. 

The  prototype  environment  that  was  developed  in  Phase  I  highlighted  the  exchange  of  model 
information  among  various  modules. 

The  prototype  that  was  developed  includes  a  module  for  automating  the  generation  of  panel 
meshes  (qtiad  panels).  Meshes  can  be  generated  from  native  AML  geometries  or  imported 
geometries.  Standard  queries  and  data  structure  interfaces  have  been  developed  which  allow  the 
generation  of  geometric  input  data  for  ZONA’s  analysis  codes. 

Detailed  studies  and  reports  can  be  fully  automated.  All  data  exchange  between  the  various 
modules  (in  different  ^sciplines)  is  fully  automated  and  tracked  through  the  CCM.  User 
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interaction  with  the  system  is  facilitated  through  a  GUI  that  allows  for  the  exploration  of  design 
specifications  and  the  investigation  and  visualization  of  results. 


A.2  Functional  Specification  Overview 

To  achieve  the  technical  objectives  of  Phase  I,  TechnoSoft  Inc.  has  divided  the  tasks  into  five 
focus  areas. 

1.  Underlying  framework  development 

The  overall  functional  specifications  of  the  underlying  framework  was  reviewed  and  developed. 
The  prototype  framework  was  developed  using  the  Adaptive  Modeling  Language  (AML)  from 
TechnoSoft  and  inherited  the  various  aspects  of  its  object-oriented  environment.  The  event- 
driven  nature  and  dependency  tracking  of  AML  is  used  to  control  and  management  of  the 
various  resources  of  and  the  data  flow  among  the  application  modules  and  the  CCM. 

2.  General-purpose  modules  development 

The  prototype  environment  employed  a  modular  architecture  developed  around  a  kernel 
implemented  in  AML  to  create  an  underlying  framework  with  resource  management  and 
dependency  tracking.  Various  modules,  both  general-purpose  and  domain-specific,  were 
developed  as  an  extension  to  the  core  system.  While  the  domain-specific  modules  in  the 
framework  focused  on  one  engineering  discipline,  the  general-purpose  modules  provide 
functionality  that  is  independent  of  any  specific  engineering  domain  knowledge. 

a.  The  core  provides  fonctionality  inherent  within  the  AML  environment  including:  1)  a  visual 
class  and  part  model  builder,  2)  a  fully-portable  graphical  user  interface  builder  supporting  a 
3D  interactive  visualizer,  3)  a  powerful  model-querying  engine,  4)  a  model  manager 
incorporating  various  utilities  for  managing  a  distributed  model  and  multiple  users,  5)  a 
model  browser  and  inspector,  and  6)  an  object-oriented,  fully-parametric,  rule-based 
geometric  modeler  with  full  access  and  querying  of  geometric  topology.  In  addition,  the 
geometric  modeler  is  fully  compliant  to  various  modeling  standards,  including  IGES  and 
proprietary  formats  such  as  SHAPES,  ACIS,  and  ParaSolid. 

b.  Feature-based  grid  generation  module.  This  prototype  module  provides  the  capability  to 
automate  the  generation  of  panel-based  mesh  models.  Controls  for  panel  size  and 
distribution,  enabled  through  various  mechanisms,  allows  for  panel  size  and  curvature 
refinements  with  a  parametric  link  to  the  geometry.  Changes  in  part  geometry  will 
automatically  trigger  the  system  to  re-compute  the  panel  model  if  required.  The  grid 
automation  supports  the  automated  mesh  generation  of  the  outer  body  and  substructure 
meshes  required  for  structural  analysis.  The  prototype  module  developed  has  limited,  case- 
specific  functionality. 

c.  External  application  linking.  This  module  provides  a  general-purpose  environment 
supporting  the  linking,  input  and  output  data  management,  and  execution  control  of 
standalone  external  applications. 

d.  Criteria  management  and  design  exploration  and  optimization  (CMDEO)  module.  This 
module  supports  the  exploration  and  evaluation  of  the  design  using  various  MDO  methods, 
design  of  experiments,  and  enumeration  methods.  Visualization  methods  are  employed  to 
view  the  resulting  trade  studies  and  correlations  as  well  as  optimum  design  configurations. 


A-4 


78 


3.  Domain  specific  analysis  modules  development 

These  modules  focus  on  the  various  design  disciplines,  integrating  the  tools  and  applications 
used  in  the  disciplines  and  capturing  the  algorithms,  processes,  and  methods  used.  These 
modules  include  ZONA’s  aerodynamics  analysis,  aerothermal,  structural  thermal  analysis, 
thermal  protection  system  design  and  configuration,  structural  analysis,  engine  sizing  and 
configuration,  weight  and  closure,  and  trajectory  analysis. 

4.  Common  computational  model  development 

The  functional  specifications  of  the  environment  were  reviewed  and  the  CCM  was  developed  to 
support  the  control  and  exchange  of  data  between  the  geometry,  mesh,  and  analysis  modules. 
The  CMDEO  module  will  provide  direct  control  of  the  CCM  data  for  design  exploration  and 
optimization.  Geometric  parameters,  meshing  details,  and  application  input  and  output  are  fully 
controlled  through  the  dependency  management  and  tracking  mechanism  employed  by  the 
framework  to  enable  the  support  of  automated  studies  and  distributed  analyses  required  for  the 
assessment  of  design  criteria  for  the  computation  of  objectives. 

5.  Environment  integration  and  testing 

The  prototype  environment  modules  were  tested  and  integrated  and  a  limited  set  of  models  was 
tested. 

TechnoSofl  has  developed  a  prototype,  case-specific  environment  that  illustrates  the  great  value 
of  developing  a  fully  functional  TAV  design,  analysis,  and  evaluation  integration  environment. 

A.3  Results 

A  protot5q)e,  case-specific  environment  was  developed  to  demonstrate  the  feasibility  of  the 
architecture.  The  environment  allows  users  to  generate  geometry  decks  that  can  be  used  to  run 
ZONA’s  analysis  codes. 

The  prototype  environment  includes  a  number  of  modules; 

1 .  Fuselage  modeling 

2.  Wing  modeling 

3.  TAV  modeling 

4.  Import  Geometry 

5.  Mesh  generation 

6.  ZONA’s  geometry  deck  generation 
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Fuselage  Modeling 

This  portion  of  the  environment  allows  users  to  rapidly  model  a  TAV  fuselage.  This  conceptual 
level  model  allows  users  of  the  system  to  rapidly  evaluate  the  viability  of  a  fuselage  design  in  an 
efficient  manor. 


Figure  A.1  Fuselage  Modeling 

Figure  A.1  is  a  sample  fuselage  that  was  defined  in  the  environment.  Users  of  the  environment 
have  full  control  over  the  shape  of  the  fuselage,  and  can  drastically  change  the  shape  of  a 
fuselage  by  varying  the  control  parameters. 
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mng  Modeling 

Similar  to  the  fuselage  modeling  concept  described  previously,  users  of  the  system  can  describe 
a  wing  geometry  that  can  be  used  for  analysis. 


Figure  A.2  Wing  geometry 
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TAV  modeling 

This  allows  the  user  to  define  a  TAV  geometry  where  the  position  of  the  wing  and  tail  are  driven 
by  the  fuselage  shape  and  size.  The  AML  TAV  model  allows  users  to  explore  “what  if’ 
scenarios  in  an  efficient  manor. 


Figure  A.3  TAV  design  module 


Import  Geometry 

The  developed  environment  allows  users  to  import  IGES,  and  STEP  geometries.  The  user  can 
then  group  the  imported  geometries,  and  mesh  the  grouped  geometry.  The  mesh  can  then  be 
used  to  generate  geometry  decks  for  ZONA’s  analysis  codes. 
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Mesh  Generation 

The  mesh  generation  process  is  one  of  the  most  labor-intensive  and  time-consuming  components 
of  the  aerodynamic  design,  analysis,  and  optimization  process.  Great  benefits  can  be  achieved 
by  building  an  environment  that  allows  users  to  parametrically  link  the  mesh  to  geometry  and  to 
save  meshing  strategies  for  future  use. 

TechnoSoft  Inc.  is  developing  a  novel  approach  that  will  allow  users  to  parametrically  link  the 
TAV  geometry  to  the  mesh. 

The  focus  will  be  on  three  areas  of  research: 

•  parametric  meshing 

•  cleanup  of  imported  geometries 

•  attribute  propagation 
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Figure  A.4  Parametric  panel  based  mesh 

The  prototype  environment  developed  has  a  limited  case-specific  parametric-automatic  meshing 
functionality.  A  set  of  methods  has  been  implemented  that  will  allow  user  to  automatically  mesh 
imported  geometries  that  need  “cleanup”. 
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ZONA  *s  geometry  deck  generation 

TAV  geometry  and  a  fuselage  grid  are  used  to  generate  geometry  decks  that  can  are  used  to  run 
zona’s  analysis  codes. 

A  set  of  classes  that  will  allow  users  to  generate  complete  ZONA  analysis  code  decks  is  under 
development.  The  functionality  these  classes  will  provide  is  critical  to  closing  the  design, 
analysis,  and  evaluation  loop. 


A.4  Conclusion 

TechnoSoft  and  ZONA  Technology  have  a  rare  and  valuable  opportimity  to  develop  a  TAV 
design,  analysis,  and  evaluation  environment  that  will  allow  users  to  explore  “what  if’  scenarios 
in  an  efficient  manor,  and  that  will  facilitate  the  sharing  of  data  between  the  different  disciplines. 
The  protot3^e  environment  is  a  proof  of  concept  that  illustrates  both  the  value  and  need  for  such 
an  environment. 
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APPENDIX  B 


HYPERSONIC/SUPERSONIC  INVISCID 
SOLUTION  FOR  BLUNTED  CONES 


Blunted  Cones  in  hypersonic  flow  require  delicate  aerodynamic  treatment.  Unlike  sharp-cone 
hypersonic  aerodynamics,  conical  flow  does  not  hold  for  a  blunted-cone  whereby  its  flowfield  is 
non-uniform.  Nevertheless,  there  exist  various  approaches  to  the  blunted  cone  problem  due  to  its 
practical  importance  in  aerospace  applications.  Most  approaches  in  the  past  were 
phenomenological  in  formulation,  notably  the  Lees  modified  Newtonian  formula  (Ref  Bl) 

Cp=Cp^5/>7^^b  (B.l) 

where  Cp^  takes  up  the  normal  shock  solution  6^  is  the  local  surface  inclination  with  respect  to 
the  freestream. 

Further  generalizations  of  (B.l)  is  the  so-called  modified  Newtonian  theory  (MNT,B1)  for  blunt- 
nose  bodies  under  angles-of-attack,  a,  reads 

Cp  =  Cp^  sin  ^gq  (B'2) 

Note  that  of  (B.  1)  is  replaced  by  5eq  where 

5eq=  sin~^  {sin(f)Cosa  -  cos^cosBsina] 

is  the  local  inclination  angle  between  the  body  slope  to  the  freestream;  (f)  is  the  azimuthal  angle 
and  6  is  the  angle  measured  from  the  local  tangent  vector  to  the  x-axis  (Figure  B.  1),  i.e. 

6  =  —  where  r^  =  y^  +  ■£' 
dx 

Although  equations  (B.l)  and  (B,2)  serve  as  accurate  and  convenient  formulas,  the  validity  of 
these  equations  are  confined  to  the  blunt-nose  region. 


z 


y 

Figure  B.l  Nomenclature  Used  For  Determination  Of  Angle 
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B.l  The  Blunted  Cone  Flow  Regions 


The  blunt-cone  problem  that  we  intend  to  solve  is  not  confined  to  the  nose  region,  however. 
With  TAV/RLV  type  configurations  in  mind,  the  blunt-cone  model  is  of  further  extent  in  its 
conical  section.  From  test  data  and  previous  CFD  solutions  (Fig  B.2),  it  can  be  seen  that  such  a 
pressure  solution  has  typically  three  distinctive  regions.  Shown  in  Fig  B.2,  region  I  is  the  nose 
section,  where  Equations  (B.l),  (B.2),  and  (B.3)  would  be  valid.  Region  III  is  the  “conical 
section”  in  which  the  surface  pressure  returns  to  that  of  the  conical  flow  and  maintains  at  a 
constant  pressure  thereupon  according  to  that  of  the  corresponding  sharp  cone  value. 

Region  II  is  the  transition  zone,  where  the  highly  compressed  flow  in  Region  I  is  gradually 
recovered  toward  the  conical  flow  solution  of  Region  III.  There  are  two  demarcations  that 
feature  this  flow  region:  the  point  of  minimum  pressure  and  that  reaching  the  conical  pressure. 

To  solve  for  the  surface  pressure  distribution  starting  fi-om  the  tangential  point  (between  the  nose 
cylinder  and  the  cone)  toward  the  conical-section  surface  going  through  the  demarcation  points, 
is  not  a  trivial  matter.  For  slender  blunted-cones,  there  exist  analytical  blast-wave  solutions  due 
to  Cheng  et  al  (Ref  B2)  and  Chemyi  (Ref  B3).  For  general  blunted-cones,  there  exists  several 
CFD  solution  methods,  including  that  of  Cleary  (Refs  B4  &  B5)  and  Hamilton  et  al  (Ref  B6), 
among  others.  The  outstanding  work  of  Dejamette  and  his  associates  (Refs  B7  and  B8)  has 
provided  a  surface  streamhne  method  with  a  flow  solution  method  for  Region  II.  Although  the 
resulting  pressure  in  AEROHEAT  (Ref  B8)  is  well  correlated  with  various  CFD  results,  it  is 
nevertheless  a  curved-fitted  (by  Chebyshev  Polynomials)  correlation  with  the  known  CFD/MOC 
solutions  (Ref  B  9,  B 1 0) 

B.2  Pressure  Solutions  of  the  Transition  Zone 


The  blast-wave  solution  of  Cheng  and  Chernyi  for  slender  cones  could  be  expressed  in  terms  of  a 
similarity  solution  format,  i.e. 


where 


4  =  tan  9^ 

K  is  Tsiems  Hypersonic  Parameter,  K  = 
y  is  the  specific  heat  ratio  of  the  gas 


X  is  the  nondimensional  distance,  x 


1  1 

r  X  ^ 

Dn 

Dn  is  the  nose  diameter. 


(B.3) 
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Figure  B.2  Three  Regions  of  a  Blunted-Cone  in  Hypersonic/Supersonic  Flow 


After  substantial  correlation  effort  with  test  data,  Krasnov  (Ref  B 11)  has  arrived  at  a  modified 
formula  after  Cheng/Chernyi’s  similarity  format,  i.e. 


where  x  =  K'“x*,  a  =  1.8,  v  = -0.2  and /r  =  1.5  (B.4) 

Equation  (B.4),  under  nominal  conditions,  does  not  take  the  ;k -effect  into  account  and  is  named 

universal-curve  formula  by  Krasnov.  It  is  this  formula  we  have  adopted  for  the  pressure 
solutions  in  region  II. 

B.3  Blunted  Cone  at  Yaw  (AoA) 

For  blunted  cones  at  angle  of  attack,  we  adopted  an  “equivalent  local  cone”  concept.  That  is,  we 
evaluate  the  pressure  Cp  according  to  a  local  cone  angle,  which  is  based  on  the  expression: 

N 

Oj^a,<p)  =  0^{O,  o)+  cosn(j)  (B.5) 

n=l 

For  simplicity,  we  merely  pick  N=1  for  all  cases  of  AoA.  Improvement  in  the  pressure  solutions 
is  expected  for  N=2, 

Next,  substituting  ^^.ftom  (B.5)  into  equation  (B.3)  yields  the  pressure  coefficients,  Cp,  for  all 
cases  of  AoA. 


B-4 


88 


B.4  Total  Cp  Solution 

The  Cp  solution  covering  regions  I,  II,  and  III^  constructed  by  a  4*-order  polynomial  fit  within 
these  regions.  Shown  in  Figure  B.2  is  a  Cp  ~  X  plot,  whereby  point  1  (Xj )  is  the  tangent  point 
of  the  cone  surface  and  the  semisphere  nose  and  point  2  (X2 )  is  selected  at  X2  =  1 .0. 

To  connect  the  Cp  solution  of  Region  I  (due  to  Equations  (B.l,  B.2))  with  that  of  Region  III 
(where  Cp  m  =  Cp  sharp  cone  ),  we  impose  the  following  connectivities  for  Cp  expressions  of 
Equations  (B.4,  B.5),  i.e., 

1)  Specifying  Cp  (Xj )  and  — ^(xj  at  X  =  Xj , 

dx 

to  connect  Region  I  solution  with  Region  II  solution 


and 


dC  d^C 

2)  Specifying  Cp  (X2)~Cp  sharp  cone  and  letting  — j^(X2)= - ^(X2)=  0, where  X2~1.0 


to  connect  Region  II  with  Region  III  solution. 


dx 


dx 
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•  Project  Objective 

•  Phase  I  Work 

•  Phase  II  Proposal 

•  Discussion 
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Total  Hypersonic  Aerothermoelastic 
_ Program  Architecture _ 

-  TAV/TPS Design  Strategy 


•  Parametric  Geometry 
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Phase  II  Technical  Objectives 

Goal: 

The  ZONA  team  proposes  to  continue  the  phase  I  development  in  a  unified 
hypersonic/supersonic  aerothermoelastic  methodology,  to  enhance  its 
geometry  high-fidelity  and  to  integrate  it  with  a  TAV/TPS  structural 
design/optimization  procedure  for  TAV/TSP  weight  minimization 
throughout  TAV’s  re-entry/maneuver  flight  phases. 

Further,  the  above  development  will  be  integrated  in  a  feature-based  design 
environment  with  parametric  control  of  models  and  data  exchange  using  the 
Adaptive  Modeling  Language  (AML). 

Objectives: 

1.  Development  of  an  expedient  conceptual/preliminary  design  tool  that 
allows  varying  of  parametric  geometry  for  rapid  assessment  of  design 
concepts. 

2.  The  TPS  sizing  is  to  meet  the  thermal  protection  requirements  of  the  main 
structures  while  being  included  with  the  main  structures  and  treated 
together  as  a  part  of  the  load-carrying  structure. 

3.  The  developed  software  can  be  utilized  by  a  design  engineer  as  an  outer 
loop  design  procedure  for  trajectory  optimization,  which  in  turn  will  lead 
to  a  minimum  fuel  and  weight  TAV  design. 
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Phase  I  Achievement 


•  Developed  ZSTREAM  for  streamline  solutions 

•  Integrated  ZONA7U+ZSTREAM+SHABP  for  hypersonic  aeroheating 
analysis  and  validation 

•  Validation  cases  include: 

-CKEM  Body 
-Blunt  Cones 
-X-34  Wing-Body 

•  Nose  Heat-rate  validation  with  MINIVER  through  two  trajectories 

•  Demonstrated  protypical  TPS  sizing  procedure 

•  X-34  aeroelastic  trim  analysis  using  ASTROS* 
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ZSTREAM  for  Stream  Line  Solution 


Validation:  Cases  Studied 


•  CKEMBody:  M  =  6.0,a  =  2° 

•  15°  Blunt  Cone:  M  =  10.6,  a  =  0°,  5°,  10° 

•  X-34  Wing-Body:  M  =  6.0,  a  =  9°,  15.22° 

•  Inviscid  Aerodynamics  are  validated  with  CFL3D 

•  Aerothermodynamics/Heat  rates  are  validated  with  CFL3D  +  LATCH 

•  Some  solutions  are  validated  with  test  data 

-  15°  Blimt  Cone:  Aerodynamics  /  NASA  TND-2969 

Aeroheating  /  NASA  TND-5450 

-  X-34  Wing-Body:  Aeroheating/ AIAA  98-0881 
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Pointed-Nose  CKEM  Body:  Aerodynamics 

M«,  =  6.0,a-2® 


Inviscid  Surface  Pressure  Distribution 


Wind-Side  Inviscid  Surface  Pressure 


Streamlines  Computed  by 
ZONA7U/ZSTREAM,  M  =  6.0 


ZONA7U  CFL3D/Euler 
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Pointed-Nose  CKEM  Body:  Aerodynamics 


15”  Blunt  Cone:  Aerodynamics 

M  =  10.6,  a  =  5* 


Inviscid  Surface  Pressure  Distribution 
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15”  Blunt  Cone:  Aerodynamics 


M  =  10.6,  a  =  10* 
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Laminar  Heat  Rate;  15”  Blunt  Cone 

M„=  10.6,  a  =  0°,  P„=  2.66  Ib/ft^  T„  =  89.971°R,  =  540°R 


13:111 1  LLmxIllSn 


Laminar  Heat  Rate:  15”  Blunt  Cone 

M„=  10.6,  a  =  10°,  P«,=  2.66  T«,  =  89.971°R,  =  540°R 


X-34  Wing-Body;  Aerodynamics  (II) 


M^=  6.0,  a  =  15.22° 


ZONA7U 


CFL3D 
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Sreamlines  Computed  by  ZSTREAM  on 
the  X-34  at  6;  g  =  9°,  a  =  15.22° 
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Aeroheating  of  X-34  (I) 

=  6.0,  a  =  9.0%  h  =  183  Kft.,  Hot  Wall,  Emissivity  =  0.8,  Turbulent 


ZONA7U+SHABP 


CFL3IHLATCH 

100  200  300  400  SOO  600  700  tOO  900  1000  1100  1300 

Cut-out**  due  to  singularity 
at  stagnation  point 


Front  View 


ZPCZ6/AeroThermoPns/tVPAFB J>rajlov200l 


sZONATECHIiinOGY 


Aeroheating  of  X-34  (II) 

M„  =  6.0,  a  =  15.22%  h  =  112  Kft.,  Hot  Wall,  Emissivity  =  0.8,  Turbulent 


ZONA7U+SHABP 


CFL3D+LATCH 


300  300  700  900  1100  1300  1300  1700 


300  SOO  700  900  1100  1300  tSOO  1700 


Cut-out**  due  to  singularity 
'at  stagnation  point 


ZONA7U+SHABP 


CFL3D+LATCH 


Front  View 


Lee-Side 
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Aeroheating  of  X-34  Nose  Stagnation 


0  200  400  600  800 

Time  (s) 
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X-34  Trim  Analysis  Using  ASTROS* 


X-34  Finite  Element  Model  X-34  Stress  Distribution  at 

_  M  =  6.0,a  =  9°,h  =  183Kft 

•  X-34  ASTROS*  FEM  Model  was  converted  from  MCS/NASTRAN 
Model  provided  by  Orbital/OSC 

•  TRIM  condition  is  at  M  =  6.0,  a  =  9°,  h  =  183  Kft. 

•  TRIM  results  for  total  weight  of  16,000  lbs. 

-  N,  =  0.97g. 

-  Trailing  Edge  Flap  =  2.05® 
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Phase  I  Concluding  Remarks  and 
Recommendation  for  Phase  II  Tasks 

Central  Methodologies  (Blocks  1-5)  required  for  aerothermodynamic  optimization  were  individually 
developed  and  validated. 

•  Hypersonic  Aerodynamics  and  Aerothermodynamic  methodology  for  Blocksl/2  are  aimed  at  replacing  the 
high-level  method  CFL3D+LATCH  and  the  low-level  method  MINIVER  by  a  mid-level  method 
ZONA7U+  SHABP  which  has  been  developed  in  Phase  L 

•  For  ZONA7U+SHABP  to  generate  one  set  of  X-34  aerodynamic/heat  rates  typically  requires  10  minutes 
on  a  550  MHZ  PC,  whereas  for  CFL3D+LATCH  it  requires  30  hours. 

•  ZSTREAMvfZS  developed  in  Phase  I  to  replace  the  QUADSTREAM  streamline  generator  of  SHABP  in 
that  the  latter,  derived  from  the  Newtonian  flow  consideration,  has  the  a  stagnation-point  singularity  in  its 
streamline  solution  and  it  does  not  depend  on  freestream  Mach  number. 

•  Validation/verification  of  ZONA’s  aerothermodynamic  method  ZONA7U+SHABP  suggests  that  further 
improvement  is  needed  in  the  following:  ZONA7U  requires  higher-fidelity  upgrade  in  order  to  cope  with 
the  high  AOA  and  the  lee-side  aerodynamics  of  SHABP  needs  to  be  replaced  by  the  AEROHEAT 
methodology  in  order  to  further  improve  the  local  heat  rate  estimates. 

•  The  TPS  weight  sizing  example  shows  that  the  designed  TPS  weight  can  be  further  reduced  if  an 
automated  optimized  scheme  can  be  developed.  A  database  of  TPS  material  in  terms  of  their  thermal  and 
mechanical  properties  must  be  fully  established  in  order  to  enhance  the  capability  of  the  optimized  scheme. 

•  The  trim  solution  of  the  X-34  in  terms  of  the  flight  loads,  input  to  the  structural  FEM  within  ASTORS*, 
will  yield  shear  loads  and  shock  loads  which  will  result  in  strength  constraint  in  the  ASTROS* 
optimization  procedure. 

•  Given  trajectory  inputs,  ZONA7U+SHABP  aeroheat  solution  at  the  nose  of  X-34  was  verified  with 
previous  solutions  obtained  by  NASA.  Total  optimization  loop  including  ASTROS*  will  be  tested  next 
using  an  X-34  example  as  a  demonstration  case. 

^ZOnUTECHT 
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Features  of  Improved  ZONAIR: 

A  High-Fidelity  Unified  Hypersonic/Supersonic 
Aerodynamic/Aer  othermody  namic  T  ool 


•  Original  ZONAIR  is  High-Order  Panel,  but  only  unified  in  supersonics/subsonics 

•  Incorporate  Unified  Hypersonic/Supersonic  Methodology  of  ZONA7U  into 
ZONAIR:  Local  Pulsating  Cone  Analogy 

•  Apply  Perturbed  Euler  Formulation  (PEF)  to  ZONAIR  to  account  for: 

-  Impinged  shock  strength  and  location 

-  Cross-flow  near-normal  shock  on  Lee-side  Flow 

•  Extend  high  AoA  flow  capability  using  L.E.  vortex  roll-up  model  to  complex 
TAV  configurations 

•  Extend  ZONAIR  to  include  two-body  aerodynamics  interference  effect  in 
hypersonic  flow 
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ZONAIR®:  Affigh-FideUty 
Supersonic/Subsonic  Panel  Method 


Comparison  of  ZONAIR  and  PANAIR  Paneling  Schemes 


Similar  to  structural  FEM 
(MSC/NASTRAN),  the  entire 
configuration  is  defined  by  “grids”. 

Only  the  starting  lines  of  the  wake 
need  to  be  defined  (via  CBAR 
elements). 

PATRAN,  FEMAP,  etc.,  can  be 
employed  directly  for  pre-  and 
post-processing. _ 


The  entire  configuration  is  first  divided 
into  several  “networks”.  Each  network 
is  further  divided  by  m  x  n  set  of  grids. 
The  location  of  the  wake  surfaces  must 
be  explicitly  defined. 

No  commercially  available  software  can 
be  used  directly  for  pre-  and  post 
processing. _ 


ZONAIR(II):  Shock  Location/Strength  Correction 
by  Perturbed  Euler  Formulation  (PEF) 


ZONAIR(in):  ZONAIR  Vortex  Roll-Up 
for  High  Angles-of-Attack  Aerodynamics 


•  Essential  elements  in  ZONAIR  for  vortex  roll-up  modeling: 


-  Vortex  sheet  emerging  from  the  wing  leading  edge  and  tip 
(free  vortex  sheet) 

-  Rolled-up  core  or  spiral  region  (fed  vortex  sheet)  fed  by  the 
leading-edge  and  tip-vortex  sheets 


Free/Fed  Vortex  Sheet  Kinematics  for  Vortex  RolMIp 
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Force  and  Moment  Coefficients  of  GAF  vs 


Capability  Comparison: 
ZONAIR  vs  Other  Aerodynamic  Codes 


Code 

Method 

Compu¬ 

tational 

Efficiency 

Streamline 

Solution 

For 

Aeroheating 

Hypersonic/ 
Supersonic/ 
Subsonic 
Mach  No. 

AIC  for 
Structural 
FEM 

Geometry 

High 

Fidelity 

High 

AOA 

2  Body 
Aero 
Inter¬ 
ference 

CFL3D 

Euler/N-S 

30hrs/ 

X-34 

Yes 

All 

No 

Yes 

Yes 

PANAIR 

Potential 

20  min/ 
X-34 

No 

Supersonic/ 

Subsonic 

No 

Yes 

No 

Yes 

ZONAIR 

Potential  + 
PEF 

20  min/ 
X-34 

Yes 

All 

Yes 

Linear- 
Order  Panel 

Yes 

Yes 

Potential  + 
PEF 

10  min/ 
X-34 

Yes 

All 

Yes 

Constant- 
Order  Panel 

Yes 

Potential  + 
Empirical 

<10  min 

Newtonian 

S.L. 

Empirical  for 
hypersonics 

No 

Low-Order 

Panel 

Yes 

iiiiw 

«  10  min 

No 

No  subsonics 

i 

. - . - 1 

No 

No 

No 

No 

DATCOM 

Analytical/ 

Empirical 

«  10  min 

No 

All  i 

No 

No 

Yes 

No 

AP98 

Analytical/ 

Empirical 

«  10  min 

No 

All 

No 

No 

Yes 

No 
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AEROHEAT  for  Accurate 
Aerothermodynamic  Analysis 

•  SHABP  aeroheating  analysis  is  based  on  empirical  equations  and  is  not 
accurate. 

•  AEROHEAT  solves  the  convective-heating  equations  using  an  axisymmetric 
analogy  that  allows  any  axisymmetric  boundary  layer  method  to  be  applied 
along  an  inviscid  surface  streamline. 

•  Laminar  and  turbulent  heating  rates  are  calculated  by  relating  the  surface  skin 
friction  to  the  momentum  thickness  Reynolds  number. 

•  the  3D  effects  of  the  AEROHEAT  methodology  is  included  through  the 
streamline  metric  coefficients  which  can  be  accurately  provided  by  the 
ZSTREAM  code. 

•  The  inviscid  aerodynamic  solutions  required  by  AEROHEAT  will  be 
computed  using  ZONAIR. 

•  The  integrated  ZONAIR  +  ZATREAM  +  AEROHEAT  computational 
procedure  will  be  validated  with  the  CFL3D  +  LATCH  results  of  the  CKE)M 
body,  15°  blunt  cone  and  the  X-34  wing-body  configuration. 

0ZONATBH 


TPS  Sizing  Optimization  Using 
Complex-Variable  Differentiation  Sensitivity 

•  TPS  sizing  will  be  automated  by  developing  an  optimization  driver  of  the  MINIVER/EXITS  code. 

•  For  a  given  heat  flux  applied  on  the  outer 
boundary,  the  objective  is  to  minimize  the 
total  weight  of  the  TPS  system  while  keeping 
the  temperature  at  each  layer  (T ■)  below  their 
respective  maximum  operational  temperature, 

T  ■ 

-^ov  n 

•  Minimize:  where is  the 

'■=*  density  of  the  F*  layer. 

Subjected  to:  T^  <  i  -  1,2. ..n 

Design  variables:  h^  >  0  i  =  1,2. ..n 

•  The  complex-variable  differentiation  can  provide  “numerically  exact”  derivatives  of  a 
complicated  function. 

-The  variable  hofz.  real  function  T(h)  is  replaced  by  A  +  iAh. 

-For small d/t:  T{h+i^)=^T{h)+i^^+...  Yields:  sh~  m  J 

•  To  incorporate  the  complex  variable  technique  into  the  MINIVER/EXITS  module  for  sensitivity 
analysis  is  straightforward  simply  by  declaring  all  variables  in  the  MINIVER/EXITS  module  as 
complex  variables. 

-The  imaginary  part  of  the  thickness  input  of  MINIVER/EXITS  represents  a  small 
incremental  thickness. 


i* 

Layer  1 

\ 

f 

f 

Layer  2 

J 

L 

\ 

( 

Layer  n 

To 

Typical  TPS  Sizing  Problem 


-The  sensitivity  is  the  imaginary  part  of  the  temperature  output  divided  by  the 


incremental  thickness. 
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Temperature  and  Aeroloads  Mapping  from 


Aerodynamic  to  Structural  Grids 


•  Displacement  /  forces  mapping  between  FEM  and  aerodynamic  grids  is  an 
existing  capability  in  ASTROS*.  Four  spline  methods  are  included: 

-  infinite  plate  spline  method 

-  Thin  plate  spline  method 

-  Beam  spline  method 

-  Rigid  body  attachment  method 

•  Temperature  mapping  from  aerodynamic  to 
FEM  surface  grids  will  be  developed  using  a 
fmite-element-based  mapping  procedure 

-  assumes  a  bilinear  temperature 
distribution  over  the  aerodynamic 
quadrilateral  panel 

-  Temperatures  are  defined  at  the 

comers  of  the  panel  and  then  mapped 
to  the  surface  of  the  FEM  model  Aerodynamic  to  Structure  Grids 


Temperature  Mapping  from 
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Automated  Parametric  Mesh  Generation 
for  ASTROS*  and  ZQNAIR 

•  Automated  parametric  mesh  generation  for  MSC/NASTRAN  models  is  an 
existing  capability  in  the  Supersonic  Hypersonic  Vehicle  Design  (SHVD)  system 
developed  by  LMCO/Technosoft. 

•  Because  of  the  similarity  between  the  MSC/NASTRAN  and  ASTROS*/ZONAIR 
bulk  data  cards.  A  mesh  generator  for  ASTROS*/ZONAIR  models  can  be 
developed  with  minor  modification  to  the  SHVD  system. 

•  In  order  to  monitor  the  progress  of  the  ASTROS*  optimization  computation,  a 
real-time  graphical  capability  will  be  developed  to  display  the  design  variables  on 
the  FEM  model  along  with  the  active  constraints  at  each  optimization  iteration. 

•  If  an  optimization  solution  cannot  be  achieved,  the  graphical  capability  will  help 
the  user  to  quickly  identify  the  source  of  the  problem  and  consequently  modify 
the  optimization  problem  statement  until  an  optimum  solution  is  obtained. 
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Test  Beds  of  the  Proposed  Design  Environment 

•  Hyper-X  Configuration 


•  A  Reusable  Military  Launch  System  (RMLS) 


Demonstration  of  the  Rapid  Design  Capability  of 
the  Proposed  Design  Environment 

•  Generate  an  ASTROS*  and  ZONAIR  model  based  on  a  set  of  initially  guessed  values  of  the 
geometric  controlling  parameters. 

•  Perform  a  ZONAIR+ZSTREAM+AEROHEAT analysis  to  establish  the  aerodynamic  and 
aeroheating  database. 

•  Obtain  a  trajectory  of  the  initial  design  using  POST/OTIS 

•  Conduct  a  TPS  sizing  over  the  entire  vehicle  based  on  the  heat  rate  time  histories  computed 
hy  POST/OTIS 

•  Perform  an  ASTROS*  optimization  computation  for  an  optimum  structural  design 

•  Establish  an  updated  weight  of  the  vehicle  and  initiate  a  new  trajectory  analysis 

•  Repeat  the  above  process  rmtil  a  converged  solution  is  achieved 

•  Obtain  a  converged  solution  by  defining  a  second  set  of  geometric  controlling  parameters 

•  Compare  the  weight  and  performance  of  the  two  solutions  for  the  sensitivity  of  the  RMLS 
performance  with  respect  to  the  geometric  controlling  parameters 

•  Graphically  document  all  intermediate  solutions  during  the  iteration  phase  of  the  above  two 
converged  solutions 

•  Establish  operational  guidelines  of  the  proposed  design  environment 


•  Conduct  a  final  adjustment  of  the  software  system  based  on  the  experience  gained  during  the 


RMLS  design 

ZPCZe/AeroThermoPrcsAVPAFB j>rcs_Nov200J 


^ZOUIKHNOIJOGY 


109 


Planned  Program  Schedule 
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